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Abstract 


The  feasibility  of  using  a  sparse  array  of  mirrors  in  a  satellite  system  is  being 
evaluated.  This  study  focuses  on  the  control  laws  necessary  for  achieving  the 
minimization  of  the  effect  of  vibration  which  happens  at  the  primary  mirror  of  the 
satellite,  by  using  the  eigenstructure  assignment  technique.  White  Gaussian  noise  was 
assumed  as  the  external  input.  AFRL  has  developed  a  79  state  model  with  9  control 
inputs  and  9  outputs,  3  sparse  array  mirrors,  and  a  base.  To  extend  the  theory,  a 
simplified  model  with  2  mirrors  and  a  base  which  has  12  states,  4  control  inputs  and  8 
outputs  was  developed.  The  system  is  not  completely  controllable  nor  observable  because 
there  are  4  rigid  body  modes. 

The  control  law  starts  from  realizing  the  states  which  are  controllable  or 
uncontrollable.  Then  only  for  the  controllable  states,  the  eigenstructure  assignment 
technique  is  applied.  The  closed- loop  eigenstructure  was  determined  by  analyzing  the 
open- loop  modes.  To  find  the  best  value  for  the  eigenstructure,  Newton's  line  search 
method  was  applied. 

The  relative  motions  of  the  primary  and  secondary  mirrors  were  described  by 
wave  reflection  dynamics.  Results  show  that  this  technique  achieves  a  reasonable  amount 
of  control  use  for  the  objective  and  relates  these  result  to  those  obtained  from  linear 
quadratic  regulator  (LQR)  technique . 
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CONTROLLING  THE  PRIMARY  MIRROR  IN  A  SPACE-BASED  TELESCOPE 
UTILIZING  AN  EIGENSTRUCTURE 
ASSIGNMENT  TECHNIQUE 


I.  Introduction 


Motivation 

The  United  States  derives  great  benefit  from  space-based  assets,  and  the 
importance  of  capability  of  the  satellite  is  continually  increasing.  In  1990,  NASA 
developed  and  launched  a  2.4  meter  diameter  primary  mirror  in  a  600Km  orbit.  The  effort 
to  create  large  satellites  and  launching  them  into  space  is  not  an  easy  one.  To  launch  large 
satellites  we  have  to  pay  a  lot  of  attention  to  drag,  GNC  and  manufacturing  issues.  These 
requirements  drive  the  satellites  to  be  small.  But  small  satellites  have  their  own 
restrictions  like  narrow  scope  angles.  So  to  achieve  the  advantages  of  large  satellites,  the 
Deployable  Optical  Telescope  (DOT)  is  being  considered  by  AFRL  as  the  next 
generation  of  large  optical  satellites.  To  reduce  weight,  AFRL  has  selected  sparse  array 
aperture  which  can  achieve  larger  diameters,  instead  of  filled  aperture.  To  compensate  for 
lack  of  image  they  added  high  tech  image  processing  for  reconstructing  the  perfect  image, 
albeit  dimmer  than  the  same -size  filled  aperture. 

Background  (  Description  of  AFRL’s  Deployable  Optical  Telescope  (DOT) ) 

The  following  description  was  based  on  Development  of  a  spare- aperture  testbed 
for  opto- mechanical  control  of  space- deployable  structures  from  AFRL(ref  [1]). 
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The  DOT  system  is  a  sparse 
aperture  finite-conjugate  imaging 
system  and  is  shown  in  Figure  1  [ref 
(13)].  The  primary  mirror  is  a  three- 
element,  spherical,  reflective  sparse 
aperture  array.  The  secondary  mirror  is 
a  monolithic  element,  and  is  also  a 
spherical  surface  of  rotation.  The 
design  is  similar  in  structure  to  a  large 
space  telescope  with  a  deployable 
sparse- aperture  primary  mirror  and 
deployable  secondary.  Each  primary 
mirror  segment  has  a  spherical  radius  of 
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Figure  1.  Sparse  array  mirror 


curvature  of  5  meters  and  a  clear  aperture  of  600  mm  in  diameter.  The  vertices  of  these 
segments  are  located  in  a  radially  symmetric  configuration,  a  distance  of  550  mm  from 
the  primary  optical  axis  at  increments  of  120  degrees.  The  circumscribed  aperture  of  the 
assembled  system  is  1.7  meters,  with  a  fill  factor  of  37.4  percent  and  an  equivalent  light¬ 
collecting  aperture  of  a  1 .04  meter  diameter  monolithic  mirror. 

The  impetus  for  the  finite  conjugate  optical  design  is  one  of  cost.  For  the  lab 
design  with  a  finite  object  distance,  a  point  source  may  be  used  for  full  aperture 
illumination  of  the  primary  mirror,  mitigating  the  need  for  a  nearly  2  meter  diameter 
collimator  in  the  pseudo- star  illuminator  system.  Additionally,  the  purely  spherical 
reflectors  require  only  three  degrees  of  alignment  correction,  reducing  the  complexity  of 
the  actuation  and  control  sensor  metrology  designs.  Thus,  the  lab  design  is  representative 
of  the  space  system,  but  unique. 

A  single  lightweight  primary  mirror  segment  is  shown  in  Figure  2.  The 
construction  consists  of  a  thin  ultra-low  expansion  (UFE)  glass  face  sheet,  supported  by 
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an  array  of  composite  rods  attached  to  a  lightweight  composite  structure  to  provide 
stiffness.  The  composite  portions  of  the  structure  are  nickel  coated  to  improve  thermal 
conduction  properties  and  reduce  the  distorting  effects  of  large  thermal  gradients. 

Each  primary  mirror  segment  is  attached  to  a  supporting  deployable  boom  or 
reaction  plate,  by  three  actuator  stacks.  The  actuator  assembly  consists  of  a  low  voltage 
piezo -ceramic  stack  mounted  atop  a  micrometer  drive  with  a  DC  servomotor.  The  piezo 
actuator  has  a  range  of  30  (im  and  a  resolution  of  1  nm,  providing  high  bandwidth  control 
capability  to  compensate  for  induced  mirror  jitter  due  to  external  disturbances,  such  as 
reaction  wheels.  These  devices  are  driven  by  a  O-lOOv  analog  signal  obtained  from  a 
linear  amplifier.  The  micrometer  has  a  10- mm  range  with  a  60- nm  resolution,  providing  a 
low  bandwidth,  coarse  positioning  capability  over  a  broad  range,  to  compensate  for 


Figure  2.  Ultra -lightweight  primary  mirror  segment 

deployment  errors  in  positioning  the  mirror.  The  resolution  of  the  coarse  actuators  is 
obtained  by  use  of  an  encoder  sensor  built  into  the  servomotor,  which  provides  feedback 
signals  to  a  separate  PID  controller.  This  controller  accepts  external  commands  from  the 
control  executive  via  an  RS-232  serial  interface.  Each  assembly  is  attached  to  the  reaction 


3 


plate  via  a  screw  clamp  on  the  micrometer.  A  diaphragm  flexure  assembly,  attached  to 
the  reaction  plate  by  three  posts,  protects  the  piezo-ceramic  from  shear  forces  induced  by 
gravitational  loading  when  the  reaction  plates  are  in  the  stowed  position. 

The  reaction  plates  are  constructed  of  lightweight  composite  materials.  These  sub- 
assemblies  are  attached  to  a  main  optics  integrating  structure  via  a  stainless  steel  tie-rod 
assembly  that  exhibits  dynamic  behavior  similar  to  a  latch  mechanism. 

Problem  Statement 

Many  dynamical  systems  are  modeled  using  Newton’s  Law  or  Lagrange’s 
equations  and  perturbation  theory.  The  result  is  a  second  order  system  of  linear  constant 
coefficient  differential  equations.  This  class  of  systems  can  be  mathematically  described 
by  the  equations  of  motions 


Mx  +  Cx  +  Kx  =  Du  (1) 

where  xe  "  and  u  e  91 are  the  state  and  control  (actuator)  vectors  respectively, 

M  is  the  n  x  n  positive  definite  symmetric  mass  matrix  ,  C  is  the  n  x  n  positive  semi- 

definite  symmetric  structural  damping  matrix,  K  is  the  n  x  n  positive  semi-definite 

dx 

symmetric  stiffness  matrix,  D  is  the  n  x  m  control  influence  matrix,  and  x  =  — . 

dt 

By  using  output  feedback,  the  control  and  measurement  equation  can  be  written 
as 


y  =  C,nX 

(2) 

y=cJ 

(3) 

u  =  -Gpy-Gvy 

(4) 
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where  ye  9E"  is  the  output(sensor)  vector,  and  Cm  is  the  in  x n  output 
measurement  matrix,  and  Gp  and  Gv  are  the  n  x  in  feedback  gains.  Equation  (2),  (3)  and 

(4)  can  be  substituted  into  Equation  (1)  and  everything  can  be  taken  to  the  left  hand  side: 
Mx  +  (C  +  D  GvCm  )x+(K  +  DGpCJx  =  0  (5) 

The  problem  considered  herein  is  how  to  select  the  control  matrices  [Gv  ,G  ],  such 

that  it  minimizes  a  cost  function.  The  system  considered  contains  both  rigid  body  and 
flexible  modes.  The  control  matrices  are  selected  based  on  a  placing  the  eigenvalues  and 
eigenvectors  (hereafter  refered  to  as  the  eigenstructure)  such  that  the  achieved 
eigenstructure  minimizes  the  cost  function. 

Methodology 

The  research  for  this  thesis  includes  developing  mathematical  models 
representative  of  a  sparse  telescope  array  and  designing  a  technique  and  computer 
algorithm  for  designing  a  feedback  control  system.  The  desired  eigenvalues  and 
eigenvector  (eigenstructure)  were  found  using  Newton’ s  line  search  method  by  analyzing 
the  open- loop  system  and  the  eigenstructure  assignment  technique  is  used  to  find  the 
proper  control  gain  matrix.  To  place  realistic  limits  on  the  control  usage,  the  maximum 
element  value  of  the  control  gain  matrix  was  constrained. 

The  results  of  the  eigenstructure  assignment  technique  were  compared  with 
results  obtained  using  a  linear  quadratic  regulator  (LQR)  approach. 


Organization 

This  thesis  is  organized  around  designing  a  control  system  for  adynamic  system 
with  rigid  body  and  flexible  modes.  Figure  3  illustrates  the  overall  process 
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This  thesis  begins  with  developing  the  equations  of  motion  (EOM)  for  the 
simplified  satellite  model  at  Chapter  2.  Chapter  3  establishes  the  techniques  for  achieving 
the  proper  control  gain  for  the  eigenstructure  assignment  technique.  Chapter  3  also 
addresses  the  technique  for  systems  that  do  not  have  enough  measurement  information  to 
place  all  the  desired  eigenvalues.  These  results  are  then  compared  to  LQR  results.  Then 
Chapter  4  provides  result  applied  to  the  theoretical  model.  Lastly,  the  technique  is  applied 
to  AFRL’s  DOT  experimental  data. 
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Input  System  Matrix 


Figure  3.  Block  diagram  of  Eigenstructure  Assignment  Process 
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II.  Developing  Equations  of  Motion 


Physical  modeling  of  two  primary  mirrors  and  a  base 

The  DOT  system  is  three  dimensional.  To  demonstrate  the  technique  and  aid  in 
analysis,  a  simplified  planar  model  representative  of  a  system  with  two  primary  mirrors 
and  a  base  was  used.  For  this  model,  as  in  the  actual  structure,  all  the  elements  are  subject 
to  small  perturbations,  and  thus  the  small  angle  and  small  disturbance  approximations 
were  used. 

Each  mirror’s  movement  has  both  horizontal  and  vertical  movement.  Because  of 
the  small  angle  approximation,  horizontal  movement  will  be  negligible  with  respect  to 
vertical  movement,  thus  only  vertical  directional  movement  and  rotation  was  considered. 
As  a  result,  each  mirror  has  2  degrees -of- freedom,  namely  the  vertical  movement  and 
rotation.  The  planar  model  is  shown  in  Figure  4.  The  full  system  has  a  total  of  6  degrees- 
of- freedom.  Although  there  are  six  DOFs,  only  internal  force  actuators  are  available,  thus 
not  all  the  DOFs  can  be  controlled.  For  this  model,  as  with  DOT,  the  rigid  body  modes 
can  not  be  controlled. 

For  system  control,  there  are  two  actuators  for  each  mirror  and  each  actuator  has 
its  own  spring  damper  and  electronic  forcing  device  such  as  the  piezo  actuators  used  on 
DOT. 


Deriving  the  Equations  of  Motion. 

Mirror  “1”  and  “2”  have  mass  mi  and  rm  and  moment  of  inertia  Ii  and  F 
respectively.  Next,  define  the  rotation  angle  of  mirror  “1”  and  “2”  as  1  and  2  with 
counter-clockwise  direction  for  positive  rotations.  In  this  expression  Ai,  A2,  A3  and  A4 
stand  for  the  actuators  that  have  a  spring,  a  viscous  damper,  and  an  electronic  forcing 


device  in  parallel  as  shown  in  Figure  4.  This  is  typical  for  space  applications,  where  the 
flexible  modes  are  controlled,  independent  from  the  rigid  body  (  orbital  motion  ). 
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Figure  4.  A  Planar  Physical  Model  of  a  Sparse  Array 


The  equations  of  motion  at  each  actuator  is: 

A  =  k  D  +c  D  +  F 


l  l  l 


(6) 


where  i  =1,2. .  .4,  Dj  represents  the  relative  distance  in  the  actuator  and  Z)  = 


dD 

dt 
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For  this  analysis,  the  elements  representing  the  base,  mirrors  and  secondary  mirror  were 
assumed  as  rigid  elements. 

From  a  static  point  of  view,  for  the  mirror  to  focus  on  the  focal  point,  the  angle 
should  be  : 


(7) 

(8) 


where  h  and  L  are  defined  as  in  Figure  4. 


(h\ 


klj 


<f»  =  arctan 

7t  d) 
a  = —  — 
4  2 


To  derive  the  EOMs,  first  consider  the  free  body  diagram  for  mirror  “1”. 

If  the  mirror  “1”  and  the  base  were  disturbed  in  the  positive  direction  then  the 
transversal  movement  will  be  as  shown  in  Figure  5. 


Figure  5.  Diagram  Showing  Coordinate  Relations 


Thus  the  distance  difference  (  Z)  )  for  the  actuator  1  and  2  will  be: 
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Dx  =  z,  —  zb  / 0 !  cosa  +  (L  +  /cosa)0fo 
D2  =  z{  -  zb +lQ]cosa  +  (L- lcosa)B b 

where  1  and  0,  are  defined  in  Figure  5. 

With  this  distance  difference  and  using  Equation  (6)  each  actuator  will  have  a 
force  of : 


Aj  —  klDl  +  cxDx  +  EJ 
An  —  k2D2  +  c2D2  +  F2 


Notice  that  the  positive  distance  difference  will  produce  a  compressive  force. 
The  resulting  Equations  of  Motion  for  mirror  1  is  (from  direct  application  of 
Newton’s  Law): 


'"Ft  =“A  ~A2 

7j0,  =/cosa(Aj  -A2) 


(11) 


Likewise  for  mirror  2,  the  distance  difference  for  actuator  3  and  4  will  be: 

D3  =  z2-  zb  -/0,  cosa  -  (L- 1 cosa)B h 
D4  =  z2  —zb  +/02c osa  ~(L  +  lcosa)B b 


and  each  actuator  will  have  a  force  of: 


A3  —  k2D2  +  c2D2  +  F2 
A4  =  k4D4  +  c4D4  +  F4 


In  a  similar  fashion,  the  Equations  of  Motion  for  mirror  2  and  the  base  will  be: 


it 


(14) 


m2z2  =  -A3  -  A4 
Hlh  Zj-,  =  A,  +  At  +  A  ,  +  A4 
/202  =  /cosa(A3  —  A4) 

/fc0fe  =  [L-/cosa](A,  -A2)  +  [L  +  /cosa](A4  -Aj) 

Deriving  an  Error  Metric 

With  this  set  of  dynamic  equations,  we  define  our  output  error  (  which  consist  of 
the  two  tilt  angles  and  the  two  wavelength  errors)  as  follows. 


Figure  6.  Errors  Resulting  from  Relative  Vertical  Movement. 

The  tilt  angle  errors  and  wavelength  errors  are  developed  separately.  The  total  error  is 
then  the  sum  of  all  the  errors.  From  Figure  6  the  distance  xi  and  yi  resulting  from  vertical 
movement  is: 


x1  =  a  sin  <|> 

(15) 

V[  =  a  cos(|) 

(16) 

where  a  =  zl-  zb 

Then  the  wavelength  error  (ERi)  and  tilt  angle  error  (Ti)  for  mirror  1  is: 
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(17) 

(18) 


ERX  =  —x1  =  -a  sin  c|) 

T  _  y i  _  acosty 

a/l2  +  h2  a/l2  +  /?2 

Likewise  the  wavelength  error  (ER2)  and  tilt  angle  error  (T2)  for  mirror  2  is: 


ER2  =  -b  sin  () 
b  cos(|) 

(19) 

Tl  ylL2+h2 

where  b=  z2-  zb 

(20) 

Figure  7.  Errors  Resulting  from  Mirror  1  Rotational  Movement 

From  Figure  7  the  distance  X2  and  y2  made  by  rotational  movement  of  mirror  1  is: 

y2  =a lL2  +  h2Q1  (21) 

x2  =  y2  tantj)  (22) 
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Then  the  wavelength  error  (ERi)  and  tilt  angle  error  (Ti)  for  mirror  1  is: 


r,  =Gj 


(23) 

(24) 


Likewise  the  wavelength  error  (ER2)  and  tilt  angle  error  (T2)  for  mirror  2  is: 


(25) 

(26) 


Figure  8.  Error  Resulting  from  Base  Rotational  Movement 

From  Figure  8,  the  distance  X3  and  y3  made  by  horizontal  movement  is: 

v3=ccos(|)  (27) 

y3  =  csin  (f>  (28) 

where  c  =  (h  +  hm)Qh 

The  wavelength  error  (ERi)  and  tilt  angle  error  (Ti)  for  mirror  1  is: 
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(29) 

(30) 


ERl  =  -jc3  =  ~c  cos  (|) 

T  -  yl  -  csin(fr 

1  a/zZTZ7 

Likewise  the  wavelength  error  (ER2)  and  tilt  angle  error  (T2)  for  mirror  2  is: 


(31) 

(32) 


ER2  =  x3  =ccos(f> 

,  _  y3  _  csin(|) 

a II2  +h 2  a/l2  +  h2 


The  total  wavelength  error  and  total  tilt  angle  error  made  by  the  vertical  and 
horizontal  movement  will  be  the  sum  of  all  the  above  errors.  In  addition  to  this,  the  tilt 
angle  error  is  changed  by  rotational  movements  of  mirror  1  and  2. 

ER j  =  -(z,  -  zb )  sin  <()  +  tan(|)A I L2  +h2Ql  -  (h  +  hm  )9ft  cos<|) 

(Zi-zjcos^)  (h  +  hjeb  sin4> 

1_  a/Z^TF  "  a/Z 

ER2  =  — (z2  —  z& ) sin  <|)  -  tan () -J L2  +/r02  +(h  +  hm'fib  cos(|) 

_  (Z2-Z,)c os4>  (h  +  hm )0 h  sin  (}) 

2  '  vz^tf  '  sHr+jr 

This  represents  the  coupling  of  errors  between  wavelength  error  and  tilt  angle  error. 


(33) 

(34) 

(35) 

(36) 
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Conventional  Linearization 


For  this  analysis,  we  require  a  linearized  model  described  in  state-space  form  as 


X  =  AX  +  BU 
Y  =CX  + DU 


(37) 


where  the  state  vector  X  is  defined  as  follows: 


X  =  [Z,  Z2Zb0j  0203  Z,  Z 2  Z b  0 j  0 2  03  ]T  =  [x  x] 


From  the  Equations  (9)  and  (12),  distance  vector  can  be  described  as: 


A 

'l 

0  - 

1  -/  cosa 

0 

L  +  /cosa 

A 

l 

0  - 

1  /  cosa 

0 

L  —  l  cosa 

A 

0 

1  - 

1  0 

-/ cosa 

-L  +  /cosa 

A. 

_o 

1  - 

1  0 

/  cosa 

-L-lc  osa 

Zi 

^2 

Zb 

0, 

02 

.0*. 


=  Dx 


(38) 


In  a  similar  fashion  I)  will  be  similar  type  of  D. 


A 

d2 

d3 

A 


D„ 


(39) 


The  actuator  vector  will  be  (  using  Equation  (10)  and  (13)  ): 
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V 

0 

0 

0 

C 

0 

0 

0" 

~F\ 

a2 

0 

k  2 

0 

0 

0 

c2 

0 

0 

~D~ 

+ 

F  2 

a3 

0 

0 

k3 

0 

0 

0 

C3 

0 

D 

F  3 

a4_ 

_o 

0 

0 

k4 

0 

0 

0 

c4 

F  4 

D 

D 


+ 


FI 
F  2 
F3 
F  4 


(40) 


Combining  the  results  we  can  build  Equations  of  Motion  from  Equation  (11)  and 
(14)  in  matrix  form  as: 


Z'i 

- 1  /  OTj 

- 1  /  ml 

0 

0 

z2 

0 

0 

- 1  /  m2 

- 1  /  m2 

"A, 

Zb 

1  lmb 

1/m, 

1/m, 

1/m, 

a2 

e. 

/cosa//. 

-/  cosa//j 

0 

0 

A3 

e2 

0 

0 

/  cosa  //2 

-/cosa/  /2 

A4 

A. 

—  (L  + 1  cosa )/ 1  b 

-  (L- /cosa)/ 7, 

(L-  /cosa)//. 

(L  +  /cosa)//, 

A, 


a4 


(41) 


The  output  equation  is  defined  in  terms  of  Equations  (33),  (34),  (35)  and  (36)  as: 

7,  =  ^  T2  ER ’  ERj 

where  Ti  and  T2  represent  the  tilt  angle  error  for  the  mirror  1  and  2  and  ERi  and  ERo 
represent  wavelength  error  for  the  mirror  1  and  2. 

In  matrix  form  this  is  represented  as: 
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Y  =  C. 


(42) 


Zi 
^2 
Zb 

0! 
02 
_e*_ 

where 


COS(|) 

0 

cost]) 

1 

0 

(h  +  hm)smty 

a/l2  +h2 

a/l2  +h2 

a/l2  +h2 

0 

cost]) 

cost]) 

0 

1 

(h  +  hm)  sin 

a/l2  +h 2 

a/l2  +  h 2 

a/l2  +/r2 

—sin  <f) 

0 

sintf) 

tant|)  a/  Tj2  +  h 1 

0 

—  (/)  +  /im)  COStf) 

0 

-sintf) 

sintf) 

0 

-tantf>VT7  +  If 

(h  +  hm  )cost|) 

(43) 

Thus  the  output  equation  is  given  as: 

y  =  [C„  0]X  (44) 

where  X  =  [z,  z,  0,  0,  0,  z,  z2  z4  0,  02  0„f 


Now  in  the  state- space  form,  the  A,  B,  C  and  D  matrix  will  be: 


zeroM  Im 


DDm  ■  Acm 


D 

zero 


4x6 


&ro  I* 6 

D„ 


12x12 


(45) 


where  7(6x6)  stands  for  identity  matrix. 
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zero6x4 

DD 

L  m  J 


(46) 


C  =  [C,„  0]  (47) 

and 

D  =  zero4x4  (48) 


For  a  given  set  of  physical  parameters  of  the  mirror,  base  and  actuators, 
generation  of  the  corresponding  state- space  model  matrices  A,  B,  C  and  D  have  been 
automated  and  the  MATLAB  code  is  contained  in  Appendix  A. 


Summary 

In  this  chapter  a  simplified  two  dimensional  DOT  system  and  the  corresponding 
EOMSs  were  developed.  The  state-space  model  has  less  measurement  outputs  and 
control  inputs  than  the  number  of  states.  Chapter  3  establishes  the  techniques  for 
achieving  the  proper  gains  for  the  eigenstructure  assignment  technique.  The 
eigenstructure  assignment  technique  and  the  constraints  will  be  discussed  on  a  case  by 
case  basis. 
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III.  Methodology 


The  previous  chapter  developed  the  simplified  planar  model.  This  chapter 
develops  the  methodology  for  control.  It  begins  with  a  development  of  a  performance 
index  for  the  system.  Eigenstructure  assignment  technique  will  be  developed  in  the  case 
of  the  full-  state  feedback  and  the  output  feedback.  Then  the  LQR  methods  will  be 
reviewed  to  compare  the  result. 

Developing  the  performance  index 

Control  system  performance  is  typically  measured  with  respect  to  a  performance 
metric.  For  the  problem  herein,  a  suitable  metric  must  be  established  to  measure  the 
closed- loop  performance.  For  a  linear  system  as  given  below: 

x(t)  =  Ax{t )  +  Bu  ( t ) 

(49) 

y(t)=Cx(f)  K  ’ 

Under  proportional  output  feedback  we  get  the  closed- loop  system  shown  in 
Figure  9. 


Figure  9.  General  Closed- loop  System 

For  the  sparse  array,  the  objective  is  to  minimize  the  effect  of  disturbances  (d) 
when  r  =  0.  This  is  referred  to  as  the  regulator  problem,  i.e.  (r  =  0  above)  the  objective  of 
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the  regulator  is  to  drive  e  — >  0  and  thus  y  — >  0  (or  x  — >  0  if  full- state  feedback  C  =  I  is 
available)  .  The  effect  of  measurement  error  (m)  will  be  not  considered  in  the  research. 
One  measure  of  how  successful  the  controller  is 


Ji  = 


(50) 


where  Q  is  a  positive  semidefinite  scalar  weighting  matrix  such  as 

Q  =  diag(ql,q2,q3,....qn)  (51) 


For  a  fixed  cost  Ji  the  bigger  we  make  cp  the  smaller  Xj  must  be,  and  thus  the 
speed  of  response  will  be  improved  by  increasing  a  particular  q. 

In  general  we  can  always  speed  up  the  all  response  by  increasing  the  gain  but  it 
will  cost  more  control  energy  to  achieve.  So  we  probably  need  to  add  a  penalty  for 
control  use.  Thus  consider 


(52) 


where  R  is  positive  definite.  A  natural  choice  would  be 

R  =  ding (rt  ,r2,r3  )  (53) 

and  if  we  pick  q  large,  then  u,  would  be  small  for  a  given  J2.  Thus  we  can  see  that  by 
varying  our  choices  in  Q  and  R,  we  can  trade  off  speed  of  response  for  control  use. 

If  there  is  specified  external  input  (  d  in  Figure  (9) ),  the  objective  of  control 
system  is  minimizing  the  effect  of  the  disturbance  (noise).  For  the  case  of  white  Gaussian 
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Noise,  every  output  variance  value  from  the  noise  can  be  calculated  through  a  Lyapunov 
Equation.  If  output  variance  is  defined  such  as  Vy  =  [vvl ,  vy2 , . ,v  ]T  then  the  objective 

function  can  be  the  sum  of  each  variance.  This  thesis  used  the  sum  of  all  output  variances 
for  the  objective  function. 


Response  of  Linear  system  to  white  noise 

The  following  is  based  on  reference  (11). 
Suppose  that  x(t)  is  the  solution  of 

x(t)  =  Ax(t )  +  Bw(t ) 
x(t0)=x  o 


(54) 


Where  w(t)  is  white  noise  with  intensity  V(t)  and  xo  is  a  stochastic  variable 
independent  of  w(t),  and  mean  of  x(to)  is  equal  to  mi  and  the  variance  matrix  Qo  is 


Qo  =  -  mo  X*0  -  mo  Y  } 


(55) 


Then  x(t)  has  mean 


(56) 


and  covariance 


Rx(ti’t2)  =  eM,'~,o)Q0eAh'2~,o)  +  rn(h’h)eM,^)BV(x)BTeAT^)dx 

Jtr> 


(57) 
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The  variance  matrix  Q(t)=Rx(t,t)  satisfies  the  differential  equation  (Lyapunov) 


Q(t)  =  AQ(t)  +  Q{t)AT  +BV  {t)BT 
where  Q(t0)=Qo 


(58) 


Now  consider  the  case  where  w(t)  is  stationary  i.e.  constant  intensity.  If  A  is 
asymptotically  stable  the  variance  matrix  will  be  as  follows  as  t  increases: 


lim  <2(?)  =  lim  Q(t)  =  Q  =  V  eM  BVBT  eA'x  dr 

t - ^  00  tc\  - > — 00 


(59) 


Further,  since  Q(t)  is  the  solution  to  the  differential  equation,  its  limit  Q  must 
also  satisfy  the  differential  equation  at  steady  state.  Specifically 


AQ  +  QA T  +BVBT  =  0 


(60) 


The  Lyapunov  equation  has  a  unique  solution  for  an  asymptotically  stable  system, 
which  means  that  the  all  of  the  eigenvalues  of  matrix  A  cannot  be  zeros  [ref  (12)]. 

For  the  problem  considered  herein,  this  means  if  w(t)  is  a  stationary,  zero-  mean, 
white  Gaussian  noise,  and  is  the  only  input  to  the  system 


x(t)  =  Ax(t )  +  Bw(t ) 


(61) 


then  the  process  x(t)  will  also  be  zero-mean  Gaussian  noise  with  variance  Q(t),  given  by 


Q(t)  =  eA(,-‘o)Q0eAT(t-to)  +  f  eA(,~l) BV (x)BT ey  ('^)d% 


(62) 
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For  the  steady- state  case,  and  stationary  w(t) 


Q  =  E{x(t)xT(t)} 


(63) 


and  is  the  solution  to 


AQ  +  QA T  +BVBt  =0 


(64) 


Further,  if  y(t )  =  Cx(t)  then  E{y(t)}=  0  and 

E{y(t)yT(t)}=E{Cx(t)xT(t)CT} 

=  CE{x(t)xr  (t)]cT 

=  CQCt 


(65) 


So  the  root  mean  square  response  of  y(t)  is  the  square  root  of  the  diagonal  of  the 
expression  CQCT  .  In  this  thesis,  the  objective  function  used  is  diagonal  sum  of  CQCT  . 


J  =  YJdicig(cQCT)  (66) 

where  Q  is  the  intensity  of  the  disturbance,  assumed  to  be  zero- mean  White  Gaussian 
Noise. 
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Stabilizable/  Detectable  System 

The  goal  of  control  is  to  stabilize  the  system  with  minimum  input  control.  To 
ensure  a  solution  exists,  we  must  first  verify  that  the  system  is  both  stabilizable  and 
detectable.  For  a  system,  if  there  exist  uncontrollable  modes  and  the  uncontrollable 
modes  are  unstable  then  there  is  no  way  to  stabilize  the  system.  Likewise,  if  there  exist 
unobservable  modes  and  the  unobservable  modes  were  totally  undetectable,  then  there  is 
no  way  to  detect  the  states.  Thus  before  a  control  system  can  be  designed,  we  need  to 
verify  whether  the  system  is  stabilizable  and  detectable  or  not. 

If  the  system  is  stabilizable  and  detectable,  then  it  is  reasonable  to  find  a 
controller  to  minimize  the  cost  function  in  Equation  (66). 

The  definition  of  stabilizable  is  “if  all  the  unstable  modes  are  controllable  then  the 
system  is  said  to  be  stabilizable”.  Clearly  a  controllable  system  is  sufficient  to  be 
stabilizable  but  the  reverse  is  not. 

The  definition  of  detectable  is  “if  all  the  unstable  modes  are  observable,  then  the 
system  is  said  to  be  detectable”.  Clearly  an  observable  system  is  sufficient  to  be 
detectable  but  the  reverse  is  not  [ref  (11)]. 

Thus  to  stabilize  the  system,  all  the  unstable  modes  must  be  observable  in  the 
performance  index.  Intuitively,  if  the  unstable  modes  are  not  included  in  the  performance 
index  then  the  optimization  will  not  consider  them  in  the  optimal  solution,  even  if  they 
are  controllable. 

As  developed  in  Chapter  2,  the  planar  model  has  both  uncontrollable  and 
unobservable  modes.  These  will  need  to  be  removed  prior  to  control  design. 
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Method  to  separate  a  system  into  the  controllable  and  uncontrollable  matrix 

In  general,  a  system  can  have  both  controllable  and  uncontrollable  modes.  But  to 
control  the  system  we  need  to  separate  a  system  into  the  controllable  and  uncontrollable 
form.  If  the  uncontrollable  modes  are  stable  then  the  system  can  be  stabilized. 

For  a  given  system: 


x  =  Ax  +  Bu 
y  -Cx  +  Du 


(67) 


Define  the  controllability  matrix  as  [ref(2)] 

mc=[b  ab  a2b  ....  A"'b\  (68) 


Using  singular  value  decomposition  Mc  can  be  decomposed  as  follows  : 


TLV  =  M  c 


(69) 


where  the  T  represents  an  orthogonal  matrix  and  the  columns  of  T  correspond  to  the 
singular  values  of  Mc  in  E  .  Thus  if  Mc  has  m  nonzero  singular  values  and  the 

remaining  m+1 . . ..n  form  a  basis  for  the  rest  of  Rn. 

Then  apply  transformation  T  such  as 


A  =  T  lAT ,  B  =  T~lB  ,  C  =  CT  . 


(70) 


so  that: 


_ 

Acc 

Acu 

_ 

Bcc 

A  = 

_  0 

Auu 

,  B  = 

_  0  _ 

and  C 


Cu  ]. 


(71) 
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With  this  new  state- space  relation  the  output  feedback  system  with  u  =  —Ky  is 


(72) 


(73) 


From  the  block  diagonal  form,  it  is  clear  that  the  closed-  loop  eigenvalue  is 
eigenvalue  of  Acc-BccKCc  and  eigenvalue  of  Auu. 

It  is  also  clear  from  Equation  (72)  that  Auu  cannot  be  affected  by  any  input.  So 
Auu  results  in  completely  uncontrollable  modes.  Thus  if  any  uncontrollable  modes  are 
unstable,  i.e.  if  real  part  of  any  eigenvalue  of  Auu  is  greater  than  zero,  then  the  system  is 
not  stablilizable. 

Likewise  Ace  can  be  called  the  completely  controllable  matrix. 

So  the  system  with  (Ace,  Bee,  Ccc)  is  a  completely  controllable  system. 

The  state  x  can  be  expressed  such  as  x  =  Tx  ,  then  x  can  be  decomposed  such  as 
.  Where  xcis  controllable  states  and  xu  is  uncontrollable  states. 

Notice  that  the  transfer  functions  Ccc(sl  -  Acc)~ 1  Bcc  =  C (si  -  A)'1  B  are 
equivalent  for  the  controllable  subsystem  and  the  original  system  So  we  can  use  Acc, 
Bcc  and  Ccc  matrix  instead  of  A,  B  and  C.  For  the  model  developed  in  Chapter  2,  this 
technique  must  be  applied  to  remove  the  uncontrollable  and  unobservable  (rigid  body 
modes)  before  proceeding  to  design  the  controller. 
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Analyzing  the  Open- loop  Eigenstructure 


In  a  given  system 


x 


0 

I 

M+ 

"  0  " 

-  Ks 

-Cd 

UJ 

Bs 

(74) 


,=cH+d«  =  [c,  o]M 

where  the  system  is  proper  and  Ks  represents  normalized  stiffness  matrix,  Cd  represents 
normalized  damping  matrix,  Bs  represents  input  direction  matrix  and  Cm  represents  the 
states  output. 

For  a  given  system,  if  all  the  springs  and  dampers  are  collocated,  then  the  possible 
open-  loop  mode  shapes  can  be  found  from  the  eigenvector  set  of  Ks  or  Cd.  If  the  open- 
loop  eigenstructure  doesn’t  contain  rigid  body  modes,  then  all  the  eigenvalues  of  Ks  are 
non- zeros.  The  eigenvectors  related  with  the  zero  eigenvalues  stand  for  rigid  body  mode 
shapes  and  the  others  with  non- zero  eigenvalues  are  the  mode  shapes  that  can  be 
controlled.  By  using  the  controllable  mode  shapes  as  the  basis,  any  closed- loop  mode 
shape  can  be  obtained,  constrained  only  by  the  max  allowable  control  use.  That  means 
the  basis  of  open- loop  mode  shapes  are  the  same  as  the  basis  of  closed- loop  mode  shapes. 


Developing  eigenstructure  assignment  [ref  (11)] 

Consider  a  system  that  has  n  states,  m  control  inputs,  and  r  measurement  outputs. 
For  a  static  output  feedback  system 

u  =  -Ky  (75) 

then  the  closed-  loop  regulator  is  given  by  : 
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x  =  (A  —  BKC  )x 


(76) 


and  the  eigenstructure  is  determined  by  (A  -  BKC ) 

As  described  before,  the  system  can  have  both  controllable  and  uncontrollable 
modes  at  a  same  time.  Eigenstructure  assignment  can  be  used  for  the  controllable  modes. 

x  =  Accx  +  Bcc  u 
y  -  Cccx 

Consider  the  system  that  has  n  controllable  states.  If  the  maximum  number  of 
inputs  or  outputs  is  greater  than  n ,  then  arbitrary  n  eigenvalues  can  be  assigned.  But  if 
the  maximum  number  of  inputs  or  outputs  is  less  than  n  ,  then  some  of  the  eigenvalues 
cannot  be  assigned.  There  are  several  ways  to  assign  the  eigenstructure.  This  technique 
will  briefly  be  reviewed  [ref(l  1)]. 


Eigenstructure  assignment  technique  with  output  feedback. 

For  a  linear  system  given  by: 

x  =  Ax  +  Bu 
y  =  Cx 
u  =  -Ky 

where  xe  Rn  ,u  e  Rm  ,y  =  Rr  the  closed- loop  regulator  is  given  by 


(78) 


x  =  [ A- BKC]x  (79) 

and  the  modes  are  determined  by  the  eigenvalues  of  [A- BKC] 
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If  the  system  is  completely  controllable  and  observable,  and  B  has  full  column  rank  and 
C  has  follow  row  rank,  then  the  following  property  can  be  proven  [ref(ll)]: 

1)  The  value  of  max(m,r)  closed- loop  eigenvalue  can  be  arbitrarily  assigned 
If  n  >  max(m,r)  then  the  eigenvalue  cannot  be  assigned. 

2)  The  shape  of  the  closed- loop  eigenvectors  associated  with  the  assigned 
eigenvalues  can  be  assigned  as  well.  A  total  of  min(m,r)  elements  of  each 
assigned  eigenvector  can  be  chosen  arbitrarily. 

3)  The  eigenvector  associated  with  the  closed- loop  eigenvalue  X,  must  lie  in  the 

subspace  spanned  by 


[Xj-A]  '  B 


(80) 


The  closed- loop  eigenvalue  problem  looks  like: 


[A-BKC]vi=\vi  (81) 

and  we  wish  to  choose  K  in  order  to  achieve  the  desired  Xt ,  v,  pairs. 

The  above  equation  can  be  rearranged  as  follows 


\kiI-A]vi  =  Bwi  (82) 

where  w,  =  -KCvj 

Once  the  w,  ’s  have  been  found,  the  gain  matrix  K  is  determined  by 

K  =  -W[CV]~1  (83) 

where  W  =[vv1,w2,w3,....wj 
and  V  =  [v1?  v2,  v3, . vr] 
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This  method  allows  the  designer  to  place  r  eigenvalues  assuming  [CV]  is  invertible. 


Eigenstructure  assignment  technique  for  full-state  feedback  case. 

For  the  full-  state  feedback  case,  C  =  I  and  this  implies  r  =  n. 

This  can  be  considered  as  3  separate  cases 
m=l,  lcmcn  and  m=n 

1)  m=l 

a)  K  has  n  parameters 

b)  All  n  eigenvalues  can  be  assigned 

c)  1  element  of  each  eigenvector  can  be  assigned  (but  since  eigenvectors 
have  arbitrary  scaling  in  this  case  it’s  useless) 

2)  m=n 

a)  K  has  n“  parameters 

b)  All  n  eigenvalues  can  be  assigned  arbitrarily 

c)  All  element  of  all  eigenvector  can  be  assigned  (complete 
eigenstructure  assignment  is  possible) 

3)  lcmcn 

a)  K  has  mx  n  c  n“  parameters. 

b)  All  n  eigenvalues  can  be  arbitrarily  assigned 

c)  A  portion  of  each  eigenvector  can  be  assigned 

For  each  case 


K  =  -W[CV]~l  (84) 

where  W  =[w1,w2,w3,....wr] 
and  V  —  [Vj , v2 , v3 ,...., v,. ] 
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If  the  desired  eigenvalues  are  distinct  ‘n’  values  in  the  assignment  process,  then  the 
eigenvectors  will  necessarily  be  linearly  independent  and  V  will  be  invertible.  If  the 
desired  eigenvalues  are  complex  eigenvalues  then  the  eigenvectors  will  also  be  complex, 
which  result  in  a  complex  gain  matrix  K.  To  avoid  complex  gain  matrices,  the  following 
technique  can  be  used.  Consider  Vj  and  v2  are  complex  and  conjugate. 


KV  =  —W 

=>K[v1R+ivu,vlR  ~ivu,v 3, ,vn]  =  -[w1R+iwu,wlR-iwu,w3, . ,wj 


(85) 


Post  multiply  both  sides  by 

[1/2  -i/2  0" 

1/2  it  2  0 
_  0  0  I_ 

which  yields: 

,vu , , v3 , . ,  v„  ]  =  ~[wlR ,  wu ,  w3, . ,  wn  ]  (  gg) 

Then  K=-WV'1,  and  the  resulting  K  matrix  is  real. 

Thus  there  is  always  a  real  matrix  K  by  simply  replacing  the  complex  pair  of 
eigenvectors  with  the  real  and  imaginary  parts  of  the  complex  eigenvector. 

As  described  before  for  the  case  2  and  3  the  gains  are  not  unique  and  so  the  extra 
freedom  can  be  used  to  place  the  components  of  the  eigenvectors. 

To  do  this,  the  desired  eigen  pair  must  satisfy  Equation  (87). 


[M  -  Ah  =  Bwt 

where  X, ,  v,  are  the  desired  eigen  pair 


(87) 


32 


Every  eigenvector  must  lie  in  the  subspace  spanned  by  [A,,  /  —  A]  '  B  .  This  puts  a 
restriction  on  the  achievable  eigenvectors. 

Now  let’ s  find  the  algorithm  to  achieve  a  vector  close  to  the  desired  eigenvector  if 
the  desired  eigenvector  vid  does  not  lie  in  the  subspace.  To  do  this  let’s  define  a 

performance  index 


1  T 

Ji=-(yi-vid)  (88) 

where  vid  =  desired  eigenvector 

v(  =  achievable  eigenvector 

Pt :  =  a  positive  definite  symmetric  matrix  whose  elements  are  chosen  to 
weight  the  difference  between  v,  and  vid  for  each  element. 

Using  a  Lagrange  multiplier  method, 

J,=  1  /  2(v,.  -  vid  )Pi  (v,.  -  vid ) + D : r  ([ A,. /  -  A]v,.  -  Bw, )  (89) 

To  find  the  minimum,  take  partial  derivatives  with  respect  to  v, ,  wi  and  o,  and  set 

to  zero. 


P  v.  +[A  7-  A]u.  -P.v.,  =0 

i  i  L  i  J  i  i  id 


dJ 

_ i_ 

dw. 

I 


=  -BTVi  =0 


[A,./-  A]v,  -  Bw ,  =  0 


(90) 
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Which  can  be  written  as 


\ 

< 

A 

o 

Nr 

W,- 

>  =  < 

0 

r 

c 

\ _ 

Av,d. 

where 


EM -A] 


-B  0 

0  Bt 

0  Yk,I  ~A\ 


(  *  denotes  hermitian  transpose  ) 


(91) 


(92) 


It  can  be  shown  that  as  long  as  the  system  is  completely  controllable,  N;  will  be 
non- singular. 


r 

> 

A 

o 

wi 

ii 

i 

0 

r 

c 

\ _ 

p,v,d. 

(93) 


If  an  eigen  pair  is  not  to  be  altered  setting  vv’(  =  0  and  vid equal  to  it’s  original 
open- loop  value  assures  that  the  associated  eigen  pair  remain  in  their  open- loop 
configurations. 

So  eigenstructure  assignment  techniques  can  provide  desirable  time  response 
characteristics  (eigenvalue  placement)  and  modal  decoupling  (eigenvector  specifications) 
for  the  nominal  plant,  but  there  are  no  guarantees  in  terms  of  stability  or  performance 
robustness.  This  is  left  to  the  designer  in  choosing  the  desired  eigenstructure. 
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State -space  model  with  limited  measurement  data. 

Suppose  we  have  A ,  which  has  an  eigenvalue  set  A  =  {A,  I  A,  ,A2,  A3  ,....,A(I} 

Then  the  characteristic  polynomial  for  matrix  A  is 


si  -  A 


=  (s  ~  \ 


s  +  a  vs 


n—  1 


\2  X.s’-A3)....(.s’-A„) 

+  Cl  s  ™  ~ . +  cin 


(94) 


If  A  represents  closed- loop  system  A  matrix  which  is  A-B*K*C,  then  the 
equation  (94)  can  be  described  as  follows: 


si  -  A 


E /(*»)*' 


g=0 


(95) 


where  kt)  represents  the  ith  row  and  jth  column  element  of  regulator  gain  matrix  K. 

If  the  desired  closed- loop  poles  are  specified  as  eigenvalue  set  A  then  the 
characteristic  polynomial  for  matrix  A  must  satisfy  the  equation  (94)  and  (95). 

By  comparing  the  coefficient  of  the  characteristic  polynomial  between  equation 
(94)  and  (95),  every  element  of  the  regulator  gain  matrix  K  can  be  found.  This  method 
can  be  used  for  the  multi- input- multi- output  case  too. 

In  single-input-single-output  case  reference  (2)  page  833-834  shows  an  easy  way 
determining  the  K  matrix.  It  should  be  also  noted  that  in  the  SISO  case  the  K  matrix  is 
unique.But  in  the  multi- input- multi-output  case,  the  K  matrix  is  not  unique.  There  are 
several  ways  to  choose  the  K  matrix. 
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As  developed  in  Chapter  2,  the  model  has  4  inputs,  4  outputs  and  12  states.  The 
model  output  information  is  only  for  the  displacement  states.  The  dynamic  system  can  be 
described  as  follows: 


+  Bu 


0 

I 

jx}+ 

■  0  " 

Ks 

-Cd 

w 

Bs 

Y  =  C{[  +  Du  =  [Cm  0 


(96) 


lx)  l*J 

where  x  represents  the  displacement  state  vector  and  x  the  velocity.  If  the  system  is 
strictly  proper,  then  the  D  matrix  can  be  removed.  Here  Ks  represents  normalized 
stiffness  matrix,  Cd  represents  normalized  damping  matrix,  Bs  represents  input  direction 
for  the  forcing  element  and  Cm  represents  output. 

The  open- loop  system  characteristic  equation  is: 


det  [A  —  A/] 


-XI 

-Ks 


I 

-  Cd  -  XI 


(97) 


The  coefficients  of  Xn  and  X"  1  (the  two  highest  polynomial  terms)  happen  only  at  the 
product  of  the  diagonal  term. 

Now  let’s  assume  that  we  can  measure  only  displacement. 

y=\C.  oirj  (98) 

With  this  open- loop  system  let’s  build  the  closed- loop  system  with  output  feedback. 

u  =  -Ky  (99) 
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Then  the  system  will  be: 


x  = 


o  i  jyi  r  o 

—  Ks  -Cd_\xj  Bs 
0  / 
-Ks-BsKC„,  -Cel 

m 

A  clH 


K[Cm  0] 


|.V  | 

|i  I 


where 


A  cl= 


0 

Ks  -  BsKC„ 


I 

Cd 


(100) 


(101) 


The  characteristic  equation  of  the  closed- loop  system  A_cl  matrix  is: 


det[A_c/  -A l]  = 


-A/ 

-Ks  - BsKCm 

m 


I 

-Cd-Xl 


(102) 


Again,  the  coefficients  of  A,"  and  A,”-1  happen  only  at  the  product  of  the  diagonal 
term.  But  in  this  product  any  control  element  of  K  matrix  cannot  affect  the  coefficient  of 
the  highest  and  the  next  highest  term  of  A, .  That  means  that  the  coefficient  of  A,'1-1  is 
fixed.  But  as  noticed  the  coefficient  of  A,”-1  for  the  closed- loop  system  is  the  same  as  the 
open- loop  system.  In  calculating  the  polynomial,  the  coefficient  of  the  A,"-1  term  of  the 
characteristic  equation  is  equal  to  the  sum  of  all  eigenvalues.  If  one  of  the  desired 
eigenvalues  is  relatively  larger  than  before,  then  the  other  eigenvalue  should  be  relatively 
smaller  and  the  sum  of  eigenvalues  should  be  the  same  as  the  open- loop  eigenvalue  sum. 
If  the  controller  is  to  stabilize  a  mode  extremely  fast,  i.e.  the  desired  eigenvalue’s  real 
value  is  extremely  large  with  negative  sign,  then  at  least  one  of  the  other  eigenvalue’s 
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real  value  can  be  positive,  i.e.  the  closed- loop  system  could  become  unstable.  So  in  the 
eigenstructure  assignment,  there  is  an  inherent  trade-off  of  the  desired  eigenvalues. 


We  need  to  measure  the  output  velocity  in  order  to  guarantee  arbitrary  eigenvalue 

."CL  0 

assignment.  The  results  shown  in  Chapter  4  used 


0  C„ 


as  the  output  measurement 


matrix  C,  to  guarantee  the  placement  of  desired  eigenvalue.  In  practice,  this  could  be 
accomplished  with  both  displacement  and  velocity  sensors  or  with  velocity  sensors  and 
digital  integrators  for  the  displacement. 


Determining  the  desired  eigenvalues 

From  a  given  open-  loop  system,  the  related  frequency  response  can  be  found  by 
using  a  singular  value  plot.  The  peak  points  represent  the  output  strength  with  the  related 
natural  frequency.  The  objective  of  this  thesis  is  to  minimi/e  the  output  variance  for  a 
random  input.  If  that  output  variance  is  zero  when  the  input  is  random  (assumed  zero- 
mean  White  Gaussian  Noise),  that  is  the  best  result.  But  it  is  impossible  to  make  the 
output  zero,  but  minimizing  the  output  variance,  the  closed- loop  system  will  maintain 
focus  better  than  the  open- loop  system. 

The  output  variance  can  be  changed  by  adjusting  the  peak  singular  values  with 
related  frequency  in  open- loop  system.  If  the  peak  singular  value  is  lowered,  the  output 
variance  can  be  lowered.  The  eigenstructure  assignment  technique  can  be  used  to  lower 
the  peak  singular  values,  (which  are  directly  related  to  the  eigenvalues) 

The  idea  behind  selecting  the  lowering  of  the  peak  singular  values  is  as  follows. 
Figure  10  shows  the  relations  between  real,  imaginary  and  absolute  value  of  eigenvalue, 
where  the  open-  loop  eigenvalue  is  Xj  =  -a  +  bi . 
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Figure  10.  Relation  between  Natural  frequency  and  Damping  frequency 

The  speed  of  stabilizing  the  system  strictly  depends  on  the  real  part  of  the 
eigenvalue  which  is  called  an  attenuation  constant.  So  minimizing  the  output  variance,  or 
lowering  the  peak  singular  value,  is  maximizing  the  eigenvalue’ s  real  part  ‘a’  which 
maintaining  the  same  natural  frequency. 

But  each  singular  value  does  not  have  the  same  co-relation  with  the  output 
variance.  The  output  variance  is  dominated  by  the  lower  frequency  eigenvalues  because 
the  modes  with  lower  natural  frequencies  decay  slower  with  respect  to  higher  natural 
frequency  modes.  So  the  modes  which  occur  at  the  lowest  frequencies  and  have  high 

singular  values  are  the  best  candidate  modes  to  lower  their  peak  singular  values. 

But  selecting  the  best  real  eigenvalue  must  satisfy  o 2  <03 2  and  the  constraint  of 

the  max  allowable  gain  value.  The  control  effort  required,  which  is  directly  related  to  the 
required  gain,  is  dependant  on  the  required  changes  to  the  eigenvalue.  Thus  these  gain 
increases  were  minimized  by  increasing  the  real  part  but  maintaining  the  frequency  at 
which  the  modes  occur  [ref(ll)]. 


Determining  the  desired  eigenvector 
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As  described  in  the  section  Eigenstructure  Assignment  technique  with  output 
feedback,  the  desired  eigenvector  associated  with  the  closed- loop  eigenvalue  Xj  must  lie 

in  the  subspace  spanned  by: 


Vi_  \  =  [XlI  -  A]  '  B 


(103) 


This  is  the  first  necessary  condition  of  the  desired  eigenvector. 

In  the  full-  state  feedback  case,  or  a  totally  controllable  system,  the  condition 
stated  above  of  the  desired  eigenvector  is  always  satisfied,  but  a  system  which  has  rigid 
body  modes  or  uncontrollable  modes  is  not. 

As  described  in  the  section  on  Analyzing  the  open- loop  Eigenstructure,  every 
closed- loop  mode  shape  which  is  controllable  lies  in  the  subspace  spanned  by  the  open- 
loop  controllable  eigenvectors  (mode  shapes).  Let  M OL  be  a  basis  of  the  set  of  open- loop 
mode  shapes.  Then  the  closed- loop  mode  shape  M CL  can  be  determined  by  a  linear 
combination  of  open- loop  mode  shape. 


M  cl  ~  3/  0L  x  {c} 

{c}={a1  a2  .  akY 

where  k  is  the  number  of  linearly  independent  open-  loop  mode  shapes. 


(104) 


Then  the  closed- loop  eigenvector  VCL  with  arbitrary  eigenvalue  must  be: 


p 
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\XxM0l  x{c}_ 

L^xmolJ 

(105) 
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So  any  closed- loop  eigenvector  related  with  eigenvalue  Xj  must  lie  in  the 
Mr 


subspace  spanned  by 


OL 

XxM 


OL. 


and  let’s  call  that  as  Vi  _  2  .  This  is  the  second  necessary 


condition  of  the  desired  eigenvector. 

So  desired  eigenvector  must  satisfy  both  necessary  conditions.  But  there  is  an 
intersecting  space  which  satisfies  the  first  and  second  necessary  condition  as  shown  in 
Figure  11.  It  is  necessary  to  use  the  intersection  space  in  finding  the  desired  eigenvector. 

The  intersection  space  must  lie  in  the  subspace  of  the  first  and  second  necessary 
condition.  If  {x}  is  a  vector  which  lies  in  the  intersection  space  of  the  two  necessary 
conditions  space,  then  the  vector  {x}  can  be  expressed  such  as: 


{x}  =  [Vi_\fo}=  [ Vi  _  2\b)  (106) 

where  {a}  and  {b}  represents  linear  combination  of  column  space  of  the  intersection  of 
matrix  Vi  _  1  and  Vi_  2. 
i.e. 


[w_lXo}-[Vi_2M={0}  (107) 

so  [v;_l  Vi _ 2}J  “  U  {()}  (108) 


Therefore  vector  |  ^  j  is  null  space  of  [V7_l  Vi  _  2] . 

So  the  intersection  space  ( V  )  of  the  two  necessary  condition  spaces  is 
span§Vi_l\a}}oY  span{\Vi _2]{b}}.  The  space  V represents  the  achievable  desired 

eigenvector  space. 


V  =  span  {[Vz_ 

=  span{\Vi  _2]{b}} 


(109) 
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Figure  11.  Necessary  Condition  Set  Relationship  Diagram 


Then  the  best  desired  eigenvector  (V  ld )  is  a  vector  which  lies  in  the  space  of  V 
and  minimizes  the  as  in  following  equation: 


Yvi=Cvid 
J  =  YJxPxYvi 


(110) 


where  P  is  any  positive  semi- definite  penalty  matrix  and  V  id  is  a  vector  lies  in  the  space 


V. 

Physically,  equation  (110)  tells  us  to  select  eigenvectors  that  contribute  to  the 
weighted  cost 

In  reality  as  the  eigenvalue  moves,  the  related  eigenvector  moves  from  open- loop 
eigenvector  to  a  vector  which  lies  in  the  space  of  achievable  desired  eigenvector  (V).  But 
as  eigen- pair  (eigenstructure)  approaches  the  desired  eigenvalue  and  the  best  desired 
eigenvector,  the  elements  of  control  gain  matrix  increases.  The  allowable  maximum 
control  gain  constraints  limit  the  extent  to  which  the  desired  eigenvectors  can  be  achieved. 
In  this  thesis,  the  proper  desired  eigenstructure  which  satisfies  allowable  maximum 
control  gain  constraints  was  found  by  using  Newton’ s  Line  Search  Method.  Typically  a 
desired  eigenvector  was  found  by  finding  the  proper  ratio  value  r,  in  Equation  (111) 
between  1  through  0. 
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v«.=v  *  (111) 

where  V  ia  is  the  ith  mode  assigned  eigenvector  as  the  desired  one,  V  io  is  the  ith  mode  open- 
loop  system  eigenvector,  and  v  id  is  f h  mode  best  desired  eigenvector. 

Every  assigned  eigenvector  was  found  by  finding  related  r\  simultaneously  using  the 
Newton’s  Line  Search  Method  again. 

For  comparison  with  the  eigenstructure  assignment  technique,  an  LQR  approach 
was  also  used,  and  will  briefly  be  discussed  below. 


Linear  Quadradic  Regulator(LQR). 

LQR  is  the  well  known  in  finding  optimal  regulator.  As  shown  in  Equation  (112), 
the  LQR  problem  is  stated  as  determine  u(t)  which  minimize  performance  index: 

7  =  |  [xT  Qx  +  u  T  Ru  )dt  (112) 

where  Q  is  positive  semi  definite  and  R  is  positive  definite  matrix. 

The  solution  which  will  determine  the  control  law  u(t)  and  hence  compensator  is 
found  from  calculus  of  variation  techniques.  Using  a  Lagrange  multiplier  to  convert  this 
constrained  nth  order  minimization  to  2nth  order  unconstrained  minimization. 

In  the  steady  state  case,  LQR  problem  turns  out  to  be  the  solution  to  the  Algebraic 
Riccati  Equation  (ARE).  The  ARE  is  solved  in  MATLAB.  This  thesis  used  LQR 
command  in  MATLAB  by  assigning  the  diagonal  value  of  matrix  Q  and  R. 


Summary 
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This  chapter  established  the  technique  for  achieving  the  proper  gains  for  the 
eigenstructure  assignment  technique  and  showed  the  simple  approach  with  LQR  method. 

In  using  the  eigenstructure  assignment  technique,  the  number  of  measurement 
outputs  was  an  important  factor  in  achieving  desired  eigenvalues.  If  there  are  not  enough 
measurement  outputs  rather  than  states,  restricted  eigenvalue  could  be  assigned.  To 
overcome  this  restriction,  a  proper  state  estimator  or  more  output  measurement  was 
needed.  With  the  assumption  of  existence  of  best  estimator,  eigenstructure  assignment 
can  be  fully  developed. 

The  next  chapter  shows  the  result  using  eigenstructure  assignment  technique.  The 
result  of  using  LQR  is  also  shown  to  compare  the  result  from  eigenstructure  assignment 
technique.  Finally,  the  AFRL  model  analysis  results  are  added  at  the  end. 
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IV.  Result  and  Analysis 


The  previous  chapter  developed  eigenstructure  assignment  and  LQR  methods. 
This  chapter  shows  the  result  of  the  2-D  planar  model  with  both  techniques  and  the 
AFRL  model  analysis  result  using  the  eigenstructure  assignment  technique. 

Open-loop  system  analysis  (2-D  planar  Model) 

Assumptions: 

1) .  There  is  no  output  measurement  error 

2) .  In  developing  the  EOMS  (Equations  of  motion),  small  angle  and  small 
displacement  approximations  were  used  to  make  the  control  system  linear 

3) .  Both  displacements  and  velocities  can  be  measured. 

Constraints : 

1).  Control  gain  max  value  of  100  (max  absolute  value  of  any  element  of  ‘K’ 

matrix). 

As  shown  in  Chapter  1  problem  statement,  the  given  system  can  be  described  as 
follows: 


Mx  +  Cx  +  Kx  =  Du  (1) 

where 

M  :  Mass  matrix 
C :  Damping  matrix 
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K  :  Stiffness  matrix 


For  the  EOM  above,  the  state- space  system  can  be  developed  (with  the 
assumption  of  output  velocity  measurements)  as  follows: 
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where  Cm  represents  measurement  matrix  for  state  the  x. 

In  this  thesis,  parameters  in  the  M  ,  K  and  C  were  chosen  to  get  similar 
characteristics  of  performance  of  the  AFRL  DOT  model  system.  The  AFRL’s  DOT 


Frequency  (rad/sec) 


Figure  12.  Open- Loop  System  Singular  Plot  (  AFRL  DOT  System) 
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model  singular  value  plot  is  shown  in  Figure  12,  and  the  planar  model  singular  value  plot 
is  shown  in  Figure  13. 


Figure  13.  Open- Loop  System  Singular  Plot  (2-D  Planar  Model) 


Table  1.  Open- Loop  System  Mode 


1st  mode 

2nd  mode 

3rd  mode 

4th  mode 

5th  mode 

6th  mode 

Eigenvalue 

(xlO3) 

-0.1401 

±  2.643 li 

-0.1040 

2.2780i 

-0.0007 

±0.1812i 

-0.0009 

±0.2165i 

0 

0 

Damping 

5.29e-002 

4.56e-002 

3.62e-003 

4.33e-003 

Mode 

Shape 

Zl 

Z2 

Zb 

0i 

02 

0b 

0.6803 

-0.6803 

0 

-0.0013 

-0.0013 

0.2729 

0.7068 

0.7068 

-0.0283 

0 

0 

0 

0 

0 

0 

0.7071 

-0.7071 

0 

0.2699 

-0.2699 

0 

0.6256 

0.6256 

-0.2674 

0.  5774 
0.5774 
0.5774 

0 

0 

0 

0.4472 

0.4472 

0 

0.4472 

0.4472 

0.4472 

Controllability 

Controllable 

Uncontrollable 
(Rigid  body  mode) 
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The  open-  loop  mode  shapes  were  determined  by  finding  the  eigenvalues  and 
eigenvectors.  Each  mode’s  eigenstructure  is  shown  in  Table  1  and  each  mode  shape  is 
shown  in  Figure  14. 


Mode  3  Mode  4 

Figure  14.  Open- Loop  System  Mode  Shape  (2-D  Planar  Model) 


As  assumed  from  the  beginning,  the  satellite  system  can  have  several  disturbance 
sources  (noise).  All  the  sources  of  noise  were  assumed  to  be  White  Gaussian  Noise  with 
intensity  2x  1CT7  to  be  close  to  the  AFRL  model  which  had  nano  range  disturbances. 
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Using  the  Lyapunov  equation,  the  steady- state  output  variance  and  Root  Mean 
Square  (RMS)  value  was  computed.  There  are  4  displacements  and  angle  errors  and  4 
velocity  errors.  The  velocity  error  variance  was  not  used  in  the  objective  function. 
Intuitively,  the  tilt  angle  errors  are  more  important  than  the  displacement  errors,  so  the 
objective  function  used  was  as  follows: 

J  =  RT  xPxR  (113) 

where 

R  :  RMS  of  {7j  T2  Dl  D2f 
Ti  :  Tilt  Angle  Error  for  Mirror  1 
T2  :  Tilt  Angle  Error  for  Mirror  2 
Di  :  Displacement  Error  for  Mirror  1 
D2  :  Displacement  Error  for  Mirror  2 
P  :  penalty  positive  semi-definite  matrix 
( in  this  thesis  P  =  diag(l04  104  1  l) 


The  open- loop  RMS  results  are  provided  in  Table  2.  The  result  of  objective 
function  was  4.1718xlCT10. 


Table  2.  Open- Loop  System  Measurement  RMS  (2-D  Planar  Model) 


Planar  Moc 

lei  (  xlCT9) 

AFRL  Moc 

X 

O 

1 

CO 

Mirror 

Tilt  Error 

143 

143.3 

A  segment 

X-tilt 

Piston  Error 

2028.9 

46.68 

A  segment 
piston 

Mirror 

“2” 

Tilt  Error 

143 

170.5 

C  segment 

X-tilt 

Piston  Error 

2028.9 

60.58 

C  segment 
piston 
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Through  analyzing  the  open- loop  system  singular  value  plot,  the  candidate  modes 
to  be  changed  were  determined  as  the  mode  which  has  largest  singular  value  and 
minimum  absolute  real  eigenvalue. 

As  shown  in  Figure  13,  there  are  two  candidate  modes  to  be  changed.  At  first,  the 
objective  function  can  be  reduced  by  selecting  only  one  mode  and  increasing  the  real  part 
of  the  eigenvalue  which  satisfies  the  constraint  that  the  maximum  control  gain  matrix 
element  must  be  less  than  allowed  value.  The  constraint  used  was  100. 


Closed-loop  system  assigned  3rd  mode  eigenvalue 

The  3ld  mode  eigenvalue  -  0.657  ±  181i  has  the  smallest  real  value,  i.e.  the  slowest 
mode.  So  the  assumption  that  the  3  mode  has  a  large  contribution  to  the  objective 
function  result  was  made.  By  assigning  a  different  eigenvalue  real  part  and  the  same 
eigenvector  with  open- loop  system  at  the  same  frequency,  the  objective  function  result 
can  be  lowered.  But  as  the  eigenvalue  real  part  is  increased,  the  control  gain  matrix 
element  value  was  increased,  as  expected.  Thus,  there  was  a  unique  eigenvalue  that 
satisfies  the  constraint  that  the  maximum  the  control  gain  matrix  element  value  was  less 
than  100  (actually  it  was  equal  to  100). 

By  assigning  the  maximum  control  gain  matrix  element  value  as  100,  the  resulting 
singular  value  plot  with  respect  to  frequency  was  given  in  Figure  15. 

The  resulting  eigenvalue  difference  is  shown  in  Table  3. 


Table  3.  3rd  Mode  Eigenvalue  Com] 

parison. 

Open  Loop 

Closed  Loop 

Eigenvalue 

-0.657  ±  181i 

-  27. 1  ±  179i 

Damping 

3.62xl0'3 

1.50x10  1 
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Closed  Loop  Sigma  Plot 


Figure  15.  Closed- Loop  System  Singular  Plot  (  3rd  mode  changed  ) 


The  result  of  objective  function  was  2.1367  x  10  10  which  was  almost  51.2  %  of 
open- loop  result.  The  resulting  jittering  performance  is  shown  in  Table  4. 


Table  4.  Closed-loop  System  Measurement  RMS  (3rd  mode  changed) 


(  xl0“9) 

Open- Loop 

Planar  Model 

Closed-  Loop 

Planar  Model 

Mirror 

Tilt  Error 

143 

103 

Piston  Error 

2028.9 

1193 

Mirror 

Tilt  Error 

143 

103 

“2” 

Piston  Error 

2028.9 

1193 
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Closed-loop  system  assigned  4th  mode  eigenvalue 

Using  the  same  procedure,  the  4th  mode  eigenvalue  ( -  0.937  ±  216i )  has  as  small 
eigenvalue  real  part  as  3ld  mode.  That  means  that  the  4th  mode  still  has  a  contribution  to 
the  objective  function  result  as  much  as  the  3  mode.  Then  by  assigning  a  different 
eigenvalue  real  part  and  the  same  eigenvector  as  the  open- loop  system  at  the  same 
frequency,  the  objective  function  result  was  lowered  to  2.1469x10“'°  which  is  51.5  %  of 
open-  loop  result. 

Figure  16  shows  the  resulting  singular  value  plot  with  respect  to  frequency  with 
the  same  max  allowable  control  gain  matrix  element  value  constraint.  The  resulting 
eigenvalue  differences  are  shown  in  Table  5. 


Table  5.  4th  Mode  Eigenvalue  Comparison. 


Open  Loop 

Closed  Loop 

Eigenvalue 

-  0.937  ±216i 

-32.5  ±214i 

Damping 

4.33xl0'3 

1.50x10' 

Table  6.  Closed-Loop  S; 

rstem  Measurement  RMS  (4th  mode  changed) 

(  xl0“9) 

Open- Loop 

Planar  Model 

Closed-  Loop 

Planar  Model 

Mirror 

66-j^  ?? 

Tilt  Error 

143 

102 

Piston  Error 

2028.9 

1673 

Mirror 

“2” 

Tilt  Error 

143 

102 

Piston  Error 

2028.9 

1673 
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Singular  Values 


Frequency  (rad/sec) 

Figure  16.  Closed- Loop  System  Singular  Plot  (  4th  mode  changed  ) 


Closed-loop  system  assigned  3rd  and  4th  mode  eigenvalue 

The  previous  section  changed  the  eigenvalue  only  for  a  single  mode  and  showed 
the  possible  eigenvalue  range  that  can  be  achieved.  So  the  assumption  that  the  objective 
function  result  could  be  lowered  by  changing  two  eigenvalues  at  the  same  time  in  the 
achievable  range  was  made.  Newton’s  Line  Search  Method  which  is  developed  in  Matlab 
as  built-in  function  with  command  name  of  ‘ fmincon ’  was  used  in  finding  the  best  two 
eigenvalue  pair  to  minimize. 
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Table  7  and  Figure  17  show  the  optimized  result  only  by  changing  the 


eigenvalues,  but  with  the  same  open- loop  system  eigenvectors. 


Table  7.  3rd  and  4th  Mode  Eigenvalue  Comparison. 


Open  Loop 

Closec 

Loop 

3rd  mode 

4th  mode 

3rd  mode 

4th  mode 

Eigenvalue 

-  0.657  ±181i 

-  0.937  ±216i 

- 13.3±  181i 

- 17.4±  216i 

Damping 

3.62xl0'3 

4.33xl0'3 

7.36xl0'2 

8.02xl0'2 

Closed  Loop  Sigma  Plot 


Frequency  (rad/sec) 


Figure  17.  Closed- Loop  System  Singular  Plot 
(  3rd  and  4th  mode  changed  at  the  same  time ) 
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Table  8.  Closed-Loop  System  Measurement  RMS  (3rd  and  4th  mode  changed) 


(  xl0“9) 

Open- Loop 

Planar  Model 

Closed-  Loop 

Planar  Model 

Mirror 

Tilt  Error 

143 

32.6 

Piston  Error 

2028.9 

457.5 

Mirror 

Tilt  Error 

143 

32.6 

“2” 

Piston  Error 

2028.9 

457.5 

After  changing  only  the  eigenvalues  for  modes  3  and  4,  the  objective  function 
result  was  lowered  to  0.2164x10“'°  which  is  5.2%  of  the  open-loop  result. 


Closed-loop  system  assigned  3rd  and  4th  mode  eigenvector 


With  the  same  eigenvalues  achieved  above,  there  can  exist  desired  eigenvectors  to 
further  minimize  the  objective  function. 

Using  the  technique  from  Chapter  3  Equation  (109)  and  (110),  the  desired 
eigenvectors  were  determined  and  the  related  mode  shape  is  shown  in  Table  9. 


Table  9.  3rd  and  4th  Mode  Optimum  Mode  Shape. 


Mode 

Shape 

3rd] 

Mode 

4th  Mode 

Open  Loop 

Desired 

Open  Loop 

Desired 

Zl 

0 

-0.7054 

0.2699 

-0.7054 

Z2 

0 

-0.7054 

-0.2699 

-0.7054 

Zb 

0 

0.0282 

0 

0.0282 

0i 

0.7071 

0.0456 

0.6256 

0.0456 

02 

-0.7071 

-0.0456 

0.6256 

-0.0456 

0b 

0 

0 

-0.2674 

0 
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Mode  3  Mode  4 

Figure  18.  Desired  Closed-Loop  Mode  Shape  (2-D  Planar  Model) 


As  a  result,  using  the  desired  eigenvector  with  the  mode  shape  as  described  in 
Table  9,  the  eigenvector  was  exactly  achieved  and  the  resulting  objective  function  was 
lowered  to  0.1244x  1CT10  which  is  2.98  %  of  open- loop  result.  But  the  maximum  gain 
matrix  element  was  increased  to  9.99  x  106 .  So  by  assigning  the  desired  eigenvector 
between  open- loop  eigenvector  to  the  best  achievable  eigenvector  as  described  in 
Equaion  (113)  in  Chapter  3,  the  maximum  gain  matrix  element  could  be  lowered  to  100 
and  the  achieved  result  of  objective  function  was  0.2164x10“'°  which  is  almost  the  same 
result  of  the  previous  section.  The  achieved  closed- loop  mode  shape  was  almost  same 
with  the  result  of  previous  section  which  was  similar  to  the  open- loop  mode  shape. 

The  results  showed  that  using  the  desired  eigenvectors  can  further  reduce  the 
objective  function,  but  at  an  increase  in  gain  and  hence  control  usage.  The  technique 
allows  the  designer  to  quickly  adjust  the  amount  of  ‘shape  control’  achieved 
(eigenvector)  for  a  given  amount  of  control  use. 
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Closed-loop  control  using  LQR 

The  advantage  of  LQR  is  the  guaranteed  stability  with  plenty  of  margins.  The  cost 
function  minimized  was  Equation  (112)  Chapter  3,  repeated  here  for  ease  of  reading. 

J  =  f  (xTQx  +  utRu )dt  (]  ]2) 

The  given  system  in  this  thesis  has  uncontrollable  rigid  body  modes.  One  of  the 
requirements  in  using  LQR  is  that  the  system  must  be  completely  controllable.  So  only 
the  controllable  modes  were  used  by  separating  the  system  into  controllable  and 
uncontrollable  subsystem  as  outlined  in  Chapter  3.  Hie  LQR  problem  can  be  converted  to 
a  problem  of  finding  the  proper  penalty  value  for  the  state  and  control  usage.  To 
minimize  the  search  space,  this  thesis  assumed  that  all  the  control  usage  has  the  same 
importance  (i.e.  all  the  actuators  are  equally  weighted).  So  every  diagonal  value  of  the 
matrix  R  was  set  to  0.1.  With  this,  the  problem  is  now  finding  the  proper  diagonal  value 
of  the  Q  matrix.  To  find  the  proper  diagonal  value,  Newton’s  Line  Search  Method  was 
used  again. 

As  a  result,  the  1st  and  2nd  mode  eigenvalues  were  almost  the  same  as  the  open- 
loop  system.  The  3ld  and  4th  mode  eigenvalues  were  changed  mostly,  but  the  values  were 
slightly  different  from  the  result  of  assigning  the  3rd  and  4th  mode  eigenvalues. 

The  resulting  singular  plot  and  the  eigenvalue  comparison  is  given  in  Ligure  16 
and  in  Table  10. 

The  resulting  objective  function  was  lowered  to  0.2169  x  10  10  which  is  5.2  %  of 
open- loop  result  (same  as  the  eigenstructure  assignment  technique). 
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Closed  Loop  Sigma  Plot 


Frequency  (rad/sec) 


Figure  19.  Closed-loop  System  Singular  Plot  (Using  LQR  ) 


Table  10.  3rd  and  4th  Mode  Eigenvalue  Comparison(  Using  LQR) 


Open  Loop 

Closet 

Loop 

3rd  mode 

4th  mode 

3rd  mode 

4th  mode 

Eigenvalue 

-  0.657  ±181i 

-  0.937  ±216i 

- 13.7  ±  1 8  li 

- 16.9  ±  216i 

Damping 

3.62xlCT3 

4.33xl0"3 

7.54  x  10~2 

7.82  xlO'2 

Table  11.  Closed-Uoop  System  Measurement  RMS  (Using  UQR) 


(xl0“9) 

Open- Loop 

Planar  Model 

Closed- Loop 

Planar  Model 

Mirror 

Tilt  Error 

143 

32.6 

Piston  Error 

2028.9 

456 

Mirror 

Tilt  Error 

143 

32.6 

“2” 

Piston  Error 

2028.9 

456 

Comparison  between  Eigenstructure  assignment  technique  and  LQR  approach 
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In  many  cases  LQR  was  used  because  LQR  guarantees  the  stability  if  it  is  possible 
to  measure  the  output  exactly.  The  number  of  penalty  variables  to  be  determined  using 
LQR  was  equal  to  the  sum  of  the  number  of  states  and  the  number  of  inputs.  In  this  thesis, 
the  simplified  model  had  8  penalty  variables  and  it  was  the  same  as  the  number  of  system 
eigenvalues.  But  the  eigenstructure  assignment  technique  used  several  eigenvalues  (or 
modes)  which  have  small  real  eigenvalue  parts  and  mostly  occur  in  the  low  frequency 
band.  The  number  of  target  eigenvalues  was  less  than  number  of  states  and  should  be  less 
than  half  because  of  the  complex  conjugate  stipulation.  In  this  thesis  there  was  2 
eigenvalues  and  their  complex  conjugates,  was  used  as  the  target  mode  eigenvalue.  The 
resulting  objective  function  value  after  using  the  eigenvalue  assignment  was  pretty  much 
equivalent  with  the  result  of  using  the  LQR  method.  So  assigning  the  eigenvalue  was  a 
computationally  faster  method  than  LQR. 

Eigenvector  assignment  was  effective  to  reduce  the  objective  function  but  most  of 
the  reduction  is  attributed  to  the  eigenvalue  assignment.  Finding  the  optimum 
eigenvectors  (optimum  mass  motion)  was  the  same  as  finding  the  linear  combination 
from  the  first  and  second  desired  eigenvector  necessary  conditions.  Achieving  the 
optimized  eigenvectors  results  in  a  requirement  of  high  control  gain  values,  such  as  an 
order  of  109  in  the  simplified  model.  After  finding  and  assigning  the  best  desired 
eigenvector  between  open- loop  system  eigenvector  and  optimized  eigenvector  by 
Newtion’s  Line  Search  Method  (Equation  (111)  in  Chapter  3),  resulting  objective 
function  was  similar  to  the  result  from  assigning  just  the  eigenvalues. 
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Closed-loop  system  for  the  AFRL  DOT  system  with  Eigenvalue  Assignment 

The  actual  DOT  system  model  which  AFRL  developed  had  9  inputs  and  9  outputs 
and  79  states.  The  corresponding  open- loop  singular  value  plot  and  optical  jitter  output 
was  shown  previously  Figure  12  and  is  tabulated  in  Table  12. 

To  reduce  the  computation  time,  this  thesis  used  a  balanced  model  realization  and 
model  reduction  to  20  states,  with  9  inputs  and  9  outputs.  The  obtained  singular  value 
plot  is  shown  in  Figure  20.  After  reducing  the  system,  a  combination  of  noise  intensities 
was  assumed  for  the  reduced  model  open- loop  system  to  be  close  to  the  original  system’ s 
jitter  performance. 

The  assumed  noise  intensity  and  resulting  open- loop  performance  is  shown  in 
Tables  13  and  14. 

The  objective  function  was  determined  as  the  sum  of  output  variance  and  the 
open- loop  system  result  was  2.3629X1CT11 . 


Table  12.  Optical  Performance  (  AFRL  DOT  System  ) 


Segment 

D.O.F. 

Open- Loop  Jitter 

Specification 

Units  (rms) 

A 

piston 

46.68 

14 

nm 

X-tilt 

168.44 

95 

nrad 

X-tilt 

143.03 

95 

nrad 

B 

piston 

45.64 

14 

nm 

X-tilt 

106.57 

95 

nrad 

X-tilt 

89.96 

95 

nrad 

C 

piston 

60.58 

14 

nm 

X-tilt 

170.5 

95 

nrad 

X-tilt 

117.88 

95 

nrad 
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In  using  the  eigenstructure  assignment  technique,  a  full  state  feedback  system  or  a 


system  which  has  the  same  number  of  output  measurements  as  the  number  of  states 
guarantees  every  desired  eigenvalue  can  be  achieved.  The  system  that  AFRL  had 
developed  has  only  9  inputs  and  9  outputs.  This  is  less  than  the  number  of  states.  This 
means  the  eigenvalue  assignment  technique  has  a  restriction  (trade-off)  as  discussed  in 
Chapter  3  (State- space  model  with  limited  measurement  data).  For  this  thesis,  the 
assumption  that  this  system  had  an  optimized  state  estimator  was  made.  In  that  case, 
eigenstructure  assignment  technique  can  be  fully  implemented. 


Table  13.  Assumed  Input  Noise  Intensity 


Noise  direction  ( Input  #  ) 

1 

2 

3 

m 

5 

6 

7 

8 

9 

Intensity 

0 

131.0617 

0 

0 

0 

0 

19.9437 

322.9485 

0 

Table  14.  Optical  Performance  (Reduced  Model  Open-loop  System) 


Segment 

D.O.F. 

Open-Loop  Jitter 
(Reduced  Model) 

Open- Loop  Jitter 
(AFRL  Model) 

Units  (rms) 

A 

piston 

1773.2 

46.68 

nm 

X-tilt 

143.7 

168.44 

nrad 

Y-tilt 

144.8 

143.03 

nrad 

B 

piston 

1136.6 

45.64 

nm 

X-tilt 

105.3 

106.57 

nrad 

Y-tilt 

112.8 

89.96 

nrad 

C 

piston 

4368.3 

60.58 

nm 

X-tilt 

172.4 

170.5 

nrad 

Y-tilt 

122.1 

117.88 

nrad 
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Using  the  full-  state  feedback  system,  the  target  eigenvalue  to  be  changed  was 
determined  to  be  the  eigenvalue  which  is  located  in  the  frequency  between  262  and  447 
(  rad  /  sec  )  as  seen  in  the  Figure  20. 

After  assigning  the  desired  eigenvalue  and  using  the  open- loop  system 
eigenvector  as  the  desired  eigenvector,  the  closed- loop  system  was  developed  and  the 
resulting  objective  function  value  was  0.287  x  10  1 1  (12  %  of  open- loop  system  result) 
with  max  allowable  gain  value  1.  The  resulting  singular  value  plot  and  closed- loop  jitter 
performance  was  shown  in  Figure  21  and  table  15. 


Singular  Values 


System:  hdel 
Frequency  (rad/sec):  447 
[Singular  Values  (dB):  46.5 


System:  hdel 

Frequency  (rad/sec):  956 
Singular  Values  (dB):  26.8 


System:  hdel 
Frequency  (rad/sec):  262 
Singular  Values  (dB):  27.5 

2U  t- . / 


Frequency  (rad/sec) 


System:  hdel 
‘  Frequency  (rad/sec):  366 
Singular  Values  (dB):  45.9 


Figure  20.  Open- Loop  System  Singular  Plot  (Reduced  Model) 
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The  proper  desired  eigenvector  to  be  assigned  was  determined  by  using  the  first 
necessary  condition  which  was  described  in  Chapter  3,  and  given  below: 


Vi_  \  =  [XiI  -  A]  '  B  (103) 

The  second  necessary  condition  wasn’t  used  because  this  portion  assumed  full- 
state  feedback. 

After  finding  the  optimum  desired  eigenvectors  with  related  desired  eigenvalue, 
the  objective  function  result  which  satisfies  the  max  allowable  control  gain  matrix 
element  value  was  almost  the  same  with  the  result  of  assigning  the  desired  eigenvalue. 


Table  15.  Optical  Performance 

(Reduced  Model  Closed- loop  System  with  Assigning  Eigenvalue) 


Segment 

D.O.F. 

Closed- Loop  Jitter 
(Reduced  Model) 

Open- Loop  Jitter 
(Reduced  Model) 

Units  (rms) 

A 

piston 

598.4 

1773.2 

nm 

X-tilt 

51.6 

143.7 

nrad 

X-tilt 

52.2 

144.8 

nrad 

B 

piston 

276 

1136.6 

nm 

X-tilt 

22.5 

105.3 

nrad 

X-tilt 

33.1 

112.8 

nrad 

C 

piston 

1557.6 

4368.3 

nm 

X-tilt 

53.9 

172.4 

nrad 

X-tilt 

51.1 

122.1 

nrad 
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Singular  Values 


Figure  21.  Closed- Loop  System  Singular  Plot 
(Reduced  Model  after  Assigning  Eigenvalue) 


Summary 

In  analyzing  the  simplified  2-D  model,  the  result  between  using  eigenstructure 
assignment  technique  and  using  LQR  method  was  almost  the  same.  Both  results  showed 
94.8%  reduction  of  objective  function.  In  analyzing  the  AFRL  model,  the  objective 
function  was  reduced  to  12%  of  the  open- loop  result. 

Assigning  the  desired  eigenvectors  (i.e.  shape  control)  took  a  great  amount  of 
control  usage,  but  can  reduce  the  effect  of  noise  further.  As  described  in  assumptions,  the 
result  is  highly  dependent  on  the  value  of  max  allowable  gain  value.  So  as  the  control 
constraints  change,  the  performance  can  be  changed. 
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V.  Conclusions  and  Recommendations 


Conclusions 

This  research  effort  examined  the  use  of  eigenstructure  assignment  to  maintain 
the  focus  for  a  sparse  flexible  array  representative  of  a  future  space  imaging  system.  It 
was  expected  that  by  increasing  the  real  part  of  the  eigenvalue  of  the  target  modes  and 
controlling  the  eigenvector  (i.e.  shape  control)  to  maintain  the  focus,  eigenstructure 
assignment  technique  would  lead  to  closed- loop  control  which  minimized  the  effect  of 
disturbance.  The  solution  technique  involved  identifying  a  desired  eigenstructure,  and 
discussed  methods  of  dealing  with  uncontrollable,  unobservable  modes,  and  the  output 
feedback  (rather  than  full- state  feedback). 

The  advantage  of  the  method  was  reduction  of  the  computation  time  for  solution. 
There  were  2  parameters  needed  to  be  found  in  assigning  the  desired  eigenvalue  which 
was  the  most  effective  part  in  minimizing  the  effect  of  disturbance  rather  than  8 
parameters  in  the  LQR  method.  In  finding  the  desired  eigenvector,  a  search  method  was 
still  needed  but  could  reduce  the  computation  time  by  using  a  basis  for  the  achievable 
eigenvector  range  space  which  satisfied  the  necessary  conditions. 

In  controlling  the  AFRL  DOT  model,  which  had  less  output  measurements  than 
states,  the  existence  of  a  state  estimator  was  essential  to  achieve  the  desired  eigenvalues. 
If  there  is  an  optimum  estimator,  the  eigenstructure  assignment  technique  was  a 
computationally  faster  approach  than  LQR. 

Recommendations 

In  analyzing  the  AFRL  model,  this  thesis  assumed  full-  state  feedback  (i.e.  there  is 
optimized  and  perfect  state  estimators).  This  is  impossible  in  a  real  system.  Future 
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research  on  eigenstructure  assignment  techniques  should  take  into  account  the  existence 
of  an  estimator,  and  determine  the  effect  on  performance. 

Not  addressed  in  this  thesis  was  the  selection  of  a  desired  eigenvector  (shape)  to 
meet  a  secondary  objective  such  as  decoupling  the  ‘tip’  and  ‘tilt’  modes.  This  approach 
could  be  used  combined  with  optical  modes  (zemikes)  to  further  enhance  control  by 
decoupling  the  optical  modes. 

For  further  studies,  the  planar  model  should  be  expanded  to  a  3-D  model. 
Additionally  flexible  body  effects  of  the  individual  mirrors  could  also  be  included  to 
qualify  these  effects. 
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Appendix  A:  Main  script  file  (Eigenstructure  Assignment  Technique) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Eigenstructure  Assignment  Technique 

%  Planar  model  control  by  changing  3rd  and  4th  mode 
%  Main  script  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%driving  the  equations  of  motion  for  2  mirrors  and  a  base 
clc;clear;close  all; 


11= 

1.5; 

12= 

1.5; 

Ib= 

5; 

ml= 

l; 

m2= 

l; 

mb= 

50; 

h= 

4; 

hm= 

0.1; 

L= 

1; 

1= 

0.1; 

k= 

2.5*10A6; 

c= 

100; 

%Initial  Set  up 
phi=atan(h/L); 
al_0=pi/4- phi/2  ; 

DM=[1  0  -1  -l*cos(al_0)  0  (L+l*cos(al_0)); 

10-1  l*cos(al_0)  0  (L-l*cos(al_0)); 

0  1-10  -l*cos(al_0)  -(L-l*cos(al_0)); 

0  1-10  l*cos(al_0)  -(L+l*cos(al_0))];%distance  matrix 

DDM=[DM  zeros(4,6); 
zeros(4,6)  DM]; 

ACM=[k  0  0  0  c  0  0  0; 

OkOOOcOO; 

OOkOOOcO; 

0  0  0k0  0  0c]; 

AB=[l/ml*[-l  -10  0]; 
l/m2*[0  0-1  -1]; 
l/mb*[l  1  1  1]; 

1/1 1  *l*cos(al_0)*  [  1  -10  0]; 
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1/I2*l*cos(al_0)*[0  0  1  -1]; 

l/Ib*[(L+l*cos(al_0))*[-l  0  0  l]+(L-l*cos(al_0))*[0 -1  1  0]]]; 


A=[zeros(6)  eye(6) 

AB*ACM*DDM]; 

B=[zeros(6,4); 

AB]; 

Cm=[[cos(phi)/sqrt(LA2+hA2),  0,  -cos(phi)/sqrt(LA2+hA2),  1,0,- 
sin(phi)*(h+hm)/sqrt(LA2+hA2)  ;  %mirror  1  tilt  angle 
0,  -cos(phi)/sqrt(LA2+hA2),  cos(phi)/sqrt(LA2+hA2),  0,  1,  - 
sin(phi)*(h+hm)/sqrt(LA2+hA2)  ];  %mirror  2  tilt  angle 
-sin(phi)  0  sin(phi)  tan(phi)*sqrt(LA2+hA2)  0  - 
(h+hm)*cos(phi)  ;  %mirror  1  wavelength  error 

0  -sin(phi)  sin(phi)  0  -tan(phi)*sqrt(LA2+hA2) 

(h+hm)*cos(phi)];  %mirror  2  wavelength  error 

C=[Cm  zeros(4,6); 
zeros(4,6)  Cm]; 

D=zeros(8,4); 

[vo  do]=eig(A);do=diag(do); 

%Controllable  system  Identification 
[ABAR,BBAR,CBAR,T,KK]  =ctrbf(A,B,C); 

Ac=ABAR(5:12,5:12); 

Bc=BBAR(5:12,:); 

Cc=CBAR(:,5:12); 

Gol=ss(Ac,Bc,Cc(l:4,:),D(l:4,:)); 

[vd,lamd]  =eig(  Ac)  ;lamd=diag(lamd) ; 
wn5=abs(lamd(5));  %Natural  Frequency  (Mode  4) 
wn7=abs(lamd(7));  %Natural  Frequency  (Mode  3) 

max_k= 1 00 

x0=[real(lamd(5));real(lamd(7))]; 

%Using  Newton's  Line  Search  Method 

option=optimset('TolX',le-12,TolFun',le-20,'TolCon',le- 

18,'MaxFunEvals',10A6,'MaxIter',10A4,'Display','iter'); 

[x,fvall,exitflag, output]  = 

fmincon(  @  fun_nock,x0,  [] ,[],[],[],[],[],  @  nonlcon_  nock, option,  Ac  ,Bc  ,Cc,max_k) 

%Overlapping  Desired  Eigenvalue 
lamd(5)=x(l)+sqrt(wn5A2-x(l)A2)*i; 
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lamd(6)=x(l)-sqrt(wn5A2-x(l)A2)*i; 

lamd(7)=x(2)+sqrt(wn7A2-x(2)A2)*i; 

lamd(8)=x(2)-sqrt(wn7A2-x(2)A2)*i; 


%Eigenstructure  Assignment  Technique 
for  jj=l:8 

N=[lamd(jj)*eye(8)-Ac,  -Be  ,zeros(8,8); 
zeros(4,12)  ,Bc'; 

eye(8)  ,zeros(8,4),  (lamd(]j)*eye(8)-Ac)']; 
VWN ( :  ,jj ) =inv(N)  *  [zeros  ( 1 2 , 1 ) ;  vd( :  ,j j )  ] ; 
end 
VWN; 


for  jj=l:2:8 

VWN  ( :  ,j  j  )=real(  VWN  ( :  ,j  j ) ) ; 

VWN  ( :  ,j  j + 1  )=imag  (VWN  ( :  ,j  j + 1 ) ) ; 
end 
VWN; 

%  required  gain 

Kbar=-VWN(9:12,:)*pinv(Cc*VWN(l:8,:)) 


Sw=diag([l  1  1  l])*2*10A-7;  %  Assumed  White  Gaussian  noise  Intensity 

%open  loop  root  mean  square(lqr) 

S_x=lyap(Ac,Bc*Sw*Bc');% Variance  of  x  matrix 
S_y=Cc(l:4,:)*S_x*Cc(l:4,:)';  %Variance  of  y  matrix 
RMS_cl=diag(sqrt(S_y));  %Close  loop  Root  mean  square  of  x 

p_angle=10A4; 

P=diag([p_angle,p_angle,l,l]);  %Penalty  matrix 
f_ol=RMS_cl'  *P  *RMS_cl*  1 0 A 1 0 

%closed  loop  root  mean  square(lqr) 

S_x=lyap(Ac-Bc*Kbar*Cc,Bc*Sw*Bc');%Variance  of  x  matrix 
S_y=Cc(l:4,:)*S_x*Cc(l:4,:)';  %Variance  of  y  matrix 
RMS_cl=diag(sqrt(S_y));  %Close  loop  Root  mean  square  of  x 

p_angle=10A4; 

P=diag([p_angle,p_angle,l,l]);  %Penalty  matrix 
f_cl=RMS_cl'*P*RMS_cl*  10A 10 

%finding  Optimum  Desired  Eigenvector 
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for  ii=  1:2:3 
lamd(ii+4) 

SUB=orth(inv(  A-  lamd(ii+4)*eye(  1 2))  *B ) ; 
rank_sub=rank(S  UB ) ; 

OL=real(vo(7:12,[l,3,5,7])); 

0L_b=lamd(ii+4)*0L; 

SUB_2=[0L;0L_b] ; 

Rank_sub2=rank(SUB_2) ; 

SUB_total=[SUB  SUB_2]; 
rank_total=rank(SUB_total) ; 

xO=[l;l;l]; 

option=optimset('TolX',le-7,'TolFun',le-7,'TolCon',le-9, 
AlaxFunEvals',10A6,'MaxIter',10A4, 'Display', 'iter'); 

[x,fval,exitflag, output]  = 

fmincon(  @  fun_vec_t,xO,  [],[],[],[],[],[],  @  nonlcon_vec, option, SUB_2,C) ; 
x_f=[x;l]; 

vd_opt( :  ,ii)=SUB_2  *x_f ; 
end 

for  kk=  1:2:3 

vd_opt(:,kk)=vd_opt(:,kk)/norm(vd_opt(:,kk));  %Optimum  Eigenvector 
end 

Trans=T*vd_opt; 

%Assigning  Proper  Ratio  between  Open- Loop  Eigenvector  and  Optimum  Eigenvector 

x0=[0.5;0.5]; 

option=optimset('TolX',le-7,'TolFun',le-7,'TolCon',le- 
9,'MaxFunEvals',10A6,'MaxIter',10A4, 'Display', 'iter'); 

[x,fval,exitflag, output]  = 

fmincon(  @  fun_vec_a,xO,  [],[],[],[],[],[],  @  nonlcon_vec_a, option,  Ac,Bc,Cc,max_k,vd, lam 
d,Trans); 

vd(:,5)=x(l)*vd(:,5)+(l-x(l))*Trans(5:12,l); 
vd( :  ,6)=conj  (vd( :  ,5)) ; 

vd(:,7)=x(2)*vd(:,7)+(l-x(2))*Trans(5:12,3); 
vd(: ,  8)=conj  (vd( :  ,7)) ; 


forjj=l:8 

N=[lamd(jj)*eye(8)-Ac,  -Be  ,zeros(8,8); 
zeros(4,12)  ,Bc'; 

eye(8)  ,zeros(8,4),  (lamd(jj)*eye(8)-Ac)']; 
V WN ( : , j j ) =in v ( N)  *  [zeros  ( 1 2 , 1 ) ;  vd( :  ,j  j )  ] ; 
end 


forjj=l:2:8 
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VWN  ( :  ,j  j  )=real(  VWN  ( :  ,j  j )) ; 
VWN(:,jj+l)=imag(VWN(:,jj+l)); 
end 

%  required  gain  matrix 

Kbar=-VWN(9:12,:)*pinv(Cc*VWN(l:8,:)); 
[vc  dc]=eig(A-B  *Kbar*C);dc=diag(dc); 


%closed  loop  root  mean  square(lqr) 

S_x=lyap(Ac-Bc*Kbar*Cc,Bc*Sw*Bc');%Variance  of  x  matrix 
S_y=Cc(l:4,:)*S_x*Cc(l:4,:)';  %Variance  of  y  matrix 
RMS_cl=diag(sqrt(S_y));  %Close  loop  Root  mean  square  of  x 


f  cl=RMS  cl'*P*RMS  cl*10A10 


Appendix  B:  Objective  Function  File  (Finding  Eigenvalue) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Eigenstructure  Assignment  Technique 
%  Planar  model  control  by  3rd  and  4th  mode  change 
%  Objective  Function  File 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


function  f=fun_nock(x,Ac,Bc,Cc,max_k) 


[vd,lamd]  =eig(  Ac)  ;lamd=diag(lamd) ; 
wn5=abs(lamd(5));  %Natural  Frequency  (Mode  4) 
wn7=abs(lamd(7));  %Natural  Frequency  (Mode  3) 


lamd(5)=x(l)+sqrt(wn5A2-x(l)A2)*i; 

lamd(6)=x(l)-sqrt(wn5A2-x(l)A2)*i; 

lamd(7)=x(2)+sqrt(wn7A2-x(2)A2)*i; 

lamd(8)=x(2)-sqrt(wn7A2-x(2)A2)*i; 


for  jj=l:8 

N=[lamd(jj)*eye(8)-Ac,  -Be  ,zeros(8,8); 
zeros(4,12)  ,Bc'; 

eye(8)  ,zeros(8,4),  (lamd(jj)*eye(8)-Ac)']; 
VWN ( :  .jj)=inv(N)  *  [zeros ( 1 2 , 1 ) ;  vd( :  ,j j ) ] ; 
end 
VWN; 


for  jj=l:2:8 

VWN  ( :  ,j  j  )=real(  VWN  ( :  ,j  j )) ; 

VWN  ( :  ,j  j + 1  )=imag(  VWN  ( :  ,j  j + 1 )) ; 
end 
VWN; 

%  required  gain 

Kbar=-VWN(9:12,:)*pinv(Cc*VWN(l:8,:)); 


E=eig(Ac-Bc*Kbar*Cc);E=real(E);I=find(E>0);s=size(I); 

CON=cond(Ac-Bc*Kbar*Cc); 
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if  CON<10A9  &  s(l)==0 

Sw=diag([l  1  1  l])*2*10A-7;  %  Assumed  White  Gaussian  noise  Intensity 

%closed  loop  root  mean  square(lqr) 

S_x=lyap(Ac-Bc*Kbar*Cc,Bc*Sw*Bc');  %Variance  of  x  matrix 
S_y=Cc(l:4,:)*S_x*Cc(l:4,:)';  %Variance  of  y  matrix 
RMS_cl=diag(sqrt(S_y));  %Close  loop  Root  mean  square  of  x 


P=diag([10A4,10A4,l,l]);  %Penalty  matrix 
f=RMS_cl'*P*RMS_cl*10A10;  %Objective  Function 
else 

f=10A20; 

end 
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Appendix  C:  Constraint  Function  File  (Finding  Eigenvalue) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Eigenstructure  Assignment  Technique 
%  Planar  model  control  by  3rd  and  4th  mode  change 
%  Constraint  Function  File 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  [c,ceq]=nonlcon_nock(x,Ac,Bc,Cc,max_k) 

[vd,lamd]  =eig(  Ac)  ;lamd=diag(lamd) ; 
wn5=abs(lamd(5));  %Natural  Frequency  (Mode  4) 
wn7=abs(lamd(7));  %Natural  Frequency  (Mode  3) 

lamd(5)=x(l)+sqrt(wn5A2-x(l)A2)*i; 

lamd(6)=x(l)-sqrt(wn5A2-x(l)A2)*i; 

lamd(7)=x(2)+sqrt(wn7A2-x(2)A2)*i; 

lamd(8)=x(2)-sqrt(wn7A2-x(2)A2)*i; 


for  jj = 1 : 8 

N=[lamd(jj)*eye(8)-Ac,  -Be  ,zeros(8,8); 
zeros(4,12)  ,Bc'; 

eye(8)  ,zeros(8,4),  (lamd(jj)*eye(8)-Ac)']; 
VWN ( :  ,jj  )=inv(N)  *  [zeros  ( 1 2 , 1 ) ;  vd( :  ,j  j )  ] ; 
end 


for  jj=l:2:8 

VWN  ( :  ,j  j  )=real(  VWN  ( :  ,j  j )) ; 

VWN  ( :  ,j  j + 1 ) =imag  (VWN  ( : , j  j + 1 ) ) ; 
end 

%  required  gain 

Kbar=-VWN(9:12,:)*pinv(Cc*VWN(l:8,:)); 


for  j  j = 1 : 8 

temp(4*  (j  j- 1 )+ 1 :4*  (j  j  - 1  )+4, 1  )=Kbar( :  ,jj ) ; 
end 

temp=abs(temp); 
c=[temp-  max_k*ones(32, 1); 
x]; 

ceq=[]; 
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Appendix  D:  Objective  Function  File  (Finding  Eigenvector) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Planar  model  control  By  Eigenstructure  Assignment 

%  Objective  Function  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  f=fun_vec_a(x,Ac,Bc,Cc,max_k,vdJamd,Trans) 

vd(:,5)=x(l)*vd(:,5)+(l-x(l))*Trans(5:12,l); 
vd(:  ,6)=conj  (vd(:  ,5)) ; 

vd(:,7)=x(2)*vd(:,7)+(l-x(2))*Trans(5:12,3); 

vd(: ,  8)=conj  (vd( :  ,7)) ; 


forjj=l:8 

N=[lamd(jj)*eye(8)-Ac,  -Be  ,zeros(8,8); 
zeros(4,12)  ,Bc'; 

eye(8)  ,zeros(8,4),  (lamd(jj)*eye(8)-Ac)']; 
V WN ( :  jj  )=inv(N)  *  [zeros ( 1 2 , 1 ) ;  vd( :  ,j  j ) ] ; 
end 
VWN; 


for  jj=l:2:8 

VWN  ( :  ,j  j  )=real(  VWN  ( :  ,j  j )) ; 

VWN  ( :  ,j  j + 1 ) =imag  (VWN  ( : , j  j + 1 ) ) ; 
end 
VWN; 

%  required  gain 

Kbar=-VWN(9:12,:)*pinv(Cc*VWN(l:8,:)); 


%closed  loop  root  mean  square(lqr) 

Sw=diag([l  1  1  l])*2*10A-7;  %  Assumed  White  Gaussian  noise  Intensity 
S_x=lyap(Ac-Bc*Kbar*Cc,Bc*Sw*Bc');%Variance  of  x  matrix 
S_y=Cc(l:4,:)*S_x*Cc(l:4,:)';  %Variance  of  y  matrix 
RMS_cl=diag(sqrt(S_y));  %Close  loop  Root  mean  square  of  x 
p_angle=10A4; 

P=diag([p_angle,p_angle,l,l]);  %Penalty  matrix 
f=RMS_cr*P*RMS_cl*10A10  ; 
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Appendix  E:  Constraint  Function  File  (Finding  Eigenvector) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Planar  model  control  By  Eigenstructure  Assignment 

%  Constraint  Function  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  [c,ceq]=nonlcon_vec_a(x,Ac,Bc,Cc,max_k,vd,lamd, Trans) 

vd(:,5)=x(l)*vd(:,5)+(l-x(l))*Trans(5:12,l); 
vd( :  ,6)=conj  (vd( :  ,5)) ; 

vd(:,7)=x(2)*vd(:,7)+(l-x(2))*Trans(5:12,3); 

vd(: ,  8)=conj  (vd( :  ,7)) ; 


forjj=l:8 

N=[lamd(jj)*eye(8)-Ac,  -Be  ,zeros(8,8); 
zeros(4,12)  ,Bc'; 

eye(8)  ,zeros(8,4),  (lamd(jj)*eye(8)-Ac)']; 
VWN ( :  ,j  j  )=inv(N)  *  [zeros  ( 1 2 , 1 ) ;  vd( :  ,j  j )  ] ; 
end 


for  jj=l:2:8 

VWN  ( :  ,j  j  )=real(  VWN  ( :  ,j  j )) ; 
VWN(:,jj+l)=imag(VWN(:,jj+l)); 
end 

%  required  gain 

Kbar=-VWN(9:12,:)*pinv(Cc*VWN(l:8,:)); 
m_k=max(max(ab  s(  Kbar)) ) ; 


c=[m_k-  max_k; 
-x; 

x-ones(2,l)]; 

ceq=[]; 
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Appendix  F:  Main  script  file  (LQR) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

% 

%  Planar  model  control  By  L  Q  R  Method 
%  Main  script  file 

%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

%driving  the  equations  of  motion  for  2  mirrors  and  a  base 

clc 

clear 


close  all 

11= 

1.5; 

12= 

1.5; 

Ib= 

5; 

ml= 

l; 

m2= 

l; 

mb= 

50; 

h= 

4; 

hm= 

0.1; 

L = 

1; 

1= 

0.1; 

k= 

2.5*10A6; 

c= 

100; 

%Initial  Set  up 
phi=atan(h/L); 
al_0=pi/4- phi/2  ; 

DM=[1  0  -1  -l*cos(al_0)  0  (L+l*cos(al_0)); 

10-1  l*cos(al_0)  0  (L-l*cos(al_0)); 

0  1-10  -l*cos(al_0)  -(L-l*cos(al_0)); 

0  1-10  l*cos(al_0)  -(L+l*cos(al_0))];%distance  matrix 

DDM=[DM  zeros(4,6); 
zeros(4,6)  DM]; 

ACM=[k  0  0  0  c  0  0  0; 

OkOOOcOO; 

OOkOOOcO; 

0  0  0k0  0  0c]; 
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AB=[l/ml*[-l  -10  0]; 
l/m2*[0  0-1  -1]; 
l/mb*[l  1  1  1]; 

1/1 1  *l*cos(al_0)*  [  1  -10  0]; 

1/I2*l*cos(al_0)*[0  0  1  -1]; 

l/Ib*[(L+l*cos(al_0))*[-l  0  0  l]+(L-l*cos(al_0))*[0 -1  1  0]]]; 


A=[zeros(6)  eye(6) 

AB*ACM*DDM]; 

B=[zeros(6,4); 

AB]; 

Cm=[[cos(phi)/sqrt(LA2+hA2),  0,  -cos(phi)/sqrt(LA2+hA2),  1,0,- 
sin(phi)*(h+hm)/sqrt(LA2+hA2)  ;  %mirror  1  tilt  angle 

0,  -cos(phi)/sqrt(LA2+hA2),  cos(phi)/sqrt(LA2+hA2),  0,  1,  - 
sin(phi)*(h+hm)/sqrt(LA2+hA2)  ];  %mirror  2  tilt  angle 
[-sin(phi)  0  sin(phi)  tan(phi)*sqrt(LA2+hA2)  0  - 
(h+hm)*cos(phi)  ;  %mirror  1  wavelength  error 

0  -sin(phi)  sin(phi)  0  -tan(phi)*sqrt(LA2+hA2) 

(h+hm)*cos(phi)]];  %mirror  2  wavelength  error 

C=[Cm  zeros(4,6); 

zeros(4,6)  Cm]; 

D=zeros(8,4); 

[vo  do]=eig(A);do=diag(do); 

%Controllable  system  Identification 

[AB  AR,BB  AR,CB  AR,T,KK]  =ctrbf(A,B,C); 

Ac=ABAR(5: 12,5: 12); 

Bc=BBAR(5:12,:); 

Cc=CBAR(:,5:12); 

max_k=100 

x=[100;100;100;100;l;.l;.01;.01] 

option=optimset('TolX',le-10,'TolFun',le-19,'TolCon',le- 

10,'MaxFunEvals',10A6,'MaxIter',10A5,'Display','iter'); 

for  kkk=l:l 
load  XXX0 
x0=(x*l)/l; 

[x,fvall,exitflag, output]  = 

fmincon(  @  fun_lqr,x0,  [] ,[],[],[],[],[],  @  nonlcon_lqr,  option,  Ac,  Be,  Cc,max_k) 
%save  XXX0  x; 
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end 


Q=diag([x])*10Al; 

R=diag([l,l,l,l]*.l); 

K=lqr(Ac,Bc,Q,R); 

Kbar=K*CcA-l; 

Sw=diag([l  1  1  l])*2*10A-7;  %  Assumed  White  Gaussian  noise  Intensity 

%closed  loop  root  mean  square(lqr) 

S_x=lyap(Ac-Bc*Kbar*Cc,Bc*Sw*Bc');  %Variance  of  x  matrix 
S_y=Cc(l:4,:)*S_x*Cc(l:4,:)';  %Variance  of  y  matrix 
RMS_cl=diag(sqrt(S_y));  %Close  loop  Root  mean  square  of  x 


P=diag([10A4,10A4,l,l]);  %Penalty  matrix 

f=RMS  cl'*P*RMS  cl*  1 0A 1 0 
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Appendix  G:  Objective  Function  File  (LQR) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Planar  model  control  By  L  Q  R  Method 
%  Objective  Function  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


function  f=fun_lqr(x, Ac ,B c ,Cc ,max_k) 


Q=diag([x])*10Al;  %  State  Penalty 

R=diag([l,l,l,l]*.l);%  Control  Penlaty 
K=lqr(Ac,Bc,Q,R); 

Kbar=K*CcA-l; 

Sw=diag([l  1  1  l])*2*10A-7;  %  Assumed  White  Gaussian  noise  Intensity 

%closed  loop  root  mean  square(lqr) 

S_x=lyap( Ac- Be* Kbar*Cc,Bc*Sw*Bc');% Variance  of  x  matrix 
S_y=Cc(l:4,:)*S_x*Cc(l:4,:)';  %Variance  of  y  matrix 
RMS_cl=diag(sqrt(S_y));  %Close  loop  Root  mean  square  of  x 


P=diag([10A4,10A4,l,l]);  %Penalty  matrix 
f=RMS_cl'*P*RMS_cl*10A10; 
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Appendix  H:  Constraint  Function  File  (LQR) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Planar  model  control  By  L  Q  R  Method 
%  Constraint  Function  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


function  [c ,ceq] =nonlcon_lqr(x, Ac ,B c  ,Cc ,max_k) 


Q=diag([x])*10Al;  %  State  Penalty 

R=diag([l,l,l,l]*.l);%  Control  Penlaty 
K=lqr(Ac,Bc,Q,R); 

Kbar=K*CcA-l; 


for  jj = 1 : 8 

temp(4*  (j  j- 1 )+ 1 :  4*  (j  j  - 1  )+4, 1  )=Kbar( :  ,jj ) ; 
end 

temp=abs(temp); 
c=[temp-  max_k*ones(32, 1); 
-x+0.1*ones(8,l)]; 

ceq=[]; 


81 


Appendix  I:  Main  script  file  (AFRL  Model) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  AFRL  Model  control  By  Eigenstructure  Assignment  Technique 
%  Main  script  file 

%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% 

clc;clear;close  all; 

%  Model  Reduction 
load  tuned_best_fit_OL 
sys=d2c(fit_sys); 
n=size(sys(:,l)); 
h=minreal(sys); 

[hb,g]=balreal(h); 

index=21:79; 

hdel=modred(hb,index,'del'); 

figure(  1 ) ;sigma(hdel,  { 200,2000 } ) ;grid  on 

[A,B,C,D]=ssdata(hdel); 


%Setting  Open  Loop  Value  as  desired  one. 

[vd  lamd]=eig(A);lamd=diag(lamd); 
[xxO,I]=sort(abs(lamd));lamdd=larnd;vdd=vd; 


%Arranging  Eigenvalue  by  increasing  order 
for  kk=  1:20 

lamd(kk)=lamdd(I(kk)) ; 
vd( :  ,kk)=vdd( :  ,I(kk)) ; 
end 

for  kk=l:9 

r(kk,  1  )=real(lamd(2*kk)); 
Wn(kk,l)=abs(lamd(2*kk)); 
end 


max_k=l  %Max  allowable  gain  matrix  element. 

%Open  loop  root  mean  square(lqr) 

load  XX_modred  %  Contains  variable  for  the  White  Gaussian  noise 

intensity  (x) 
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Sw=diag(x); 


%  White  Gaussian  noise  intensity 


S_x=lyap(A,B*Sw*B'); 

S_y=C*S_x*C'; 

RMS  l=diag(sqrt(S_y)) 
f=RMS  l'*RMS  1/10A4 


%  Variance  of  x  matrix 
%  Variance  of  y  matrix 

%Closed  loop  Root  mean  square  of  x 
%Objective  Function 


%Finding  optimum  eigenvalue  assignment  with  open  loop  eigenvector  for  desired  one. 
xO=[r] 

option=optimset('TolX',le-9,'TolF  un',le-9,'TolCon',le-9, 
'MaxFunEvals',10A6,'MaxIter',40A3,'Display','iter'); 

[x,fvall,exitflag,output]  = 

fmincon(  @  fun_e_val,xO,  [] ,[] ,[] ,[]  ,[] ,[] ,  @  nonlcon_e_val,  option,  A,  B  ,C,max_k,lamd,vd, 
Wn); 

save  result_e_val  x 
for  kk=l:9 

lamd(2*kk-l)=x(kk,l)+sqrt(Wn(kk)A2-x(kk,l)A2)*i; 
lamd(2*kk)=conj(lamd(2*kk- 1)); 
end 

%  Checking  the  result  of  assigning  the  eigenvalue. 

P=eye(20);  %eigenvector  penalty  matrix 


for  j  j = 1 : 20 

N=[lamd(jj)*eye(20)-A,  -B  ,zeros(20); 
zeros(9,20+9)  ,B'; 

P  ,zeros(20,9),  (lamd(jj)*eye(20)-A)']; 
VWN(:,jj)=inv(N)*[zeros(29,l);P*vd(:,jj)]; 
end 


for  jj=l:2:20 

VWN  ( :  ,jj  )=real(  VWN  ( :  ,jj )) ; 

VWN  ( :  ,jj + 1  )=imag(  VWN  ( :  ,jj + 1 )) ; 
end 


Kbar=-VWN(20+l:20+9,:)*pinv(eye(20)*VWN(l:20,:));  %  required  gain  matrix  with 
full  state- feedback 

G_cl=ss(A-B*Kbar,B,C,D); 

figure(2)  ;sigma(G_cl,  { 200,2000 } )  ;grid  on 

%Closed  loop  root  mean  square(lqr) 

load  XX_modred  %  Contains  variable  for  the  White  Gaussian  noise 

intensity  (x) 
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Sw=diag(x);  %  White  Gaussian  noise  intensity 

S_x=lyap(A-B*Kbar,B*Sw*B');  %Variance  of  x  matrix 

S_y=C*S_x*C';  %Variance  of  y  matrix 

RMS=diag(sqrt(S_y))  %Closed  loop  Root  mean  square  of  x 

f=RMS'*RMS/10Al  1 

MAX_K=max(max(abs(Kbar)))  %Objective  Function 

save  E_YAL_result 
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Appendix  J:  Objective  Function  File(AFRL  Model) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  AFRL  Model  control  By  Eigenstructure  Assignment  Technique 
%  Objective  Function  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  f=fun_e_val(x,A,B,C,max_k,lamd,vd,Wn) 
for  kk=l:9 

lamd(2*kk-l)=x(kk,l)+sqrt(Wn(kk)A2-x(kk,l)A2)*i; 
lamd(2*kk)=conj(lamd(2*kk- 1)); 
end 

%  Checking  the  result  of  assigning  the  eigenvalue. 

P=eye(20);  %eigenvector  penalty  matrix 


for  j  j = 1 : 20 

N=[lamd(jj)*eye(20)-A,  -B  ,zeros(20); 
zeros(9,20+9)  ,B'; 

P  ,zeros(20,9),  (lamd(jj)*eye(20)-A)']; 
VWN(:,jj)=inv(N)*[zeros(29,l);P*vd(:,jj)]; 
end 


for  jj=l:2:20 

VWN  ( :  ,j  j  )=real(  VWN  ( :  ,jj)); 

VWN  ( :  ,jj + 1  )=imag(  VWN  ( :  ,jj + 1 )) ; 
end 

Kbar=-VWN(20+l:20+9,:)*pinv(eye(20)*VWN(l:20,:));  %  required  gain 


E=eig(  A-  B  *  Kbar)  ;E=real(E)  ;I=find(E>0) ;  s= size(I) ; 

C ON =cond(  A-  B  *  Kbar) ; 
if  CON<10A9  &  s(l)==0 

load  XX_modred  %  Contains  variable  for  the  White  Gaussian  noise  intensity 
Sw=diag(x);  %  White  Gaussian  noise  intensity 

%Closed  loop  root  mean  square(lqr) 

S_x=lyap(A-B*Kbar,B*Sw*B');  %Variance  of  x  matrix 
S_y=C*S_x*C';  %Variance  of  y  matrix 

RMS=diag(sqrt(S_y));  %Closed  loop  Root  mean  square  of  x 

f=RMS'*RMS/10A4; 
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else 

f=10A20; 

end 
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Appendix  K:  Constraint  Function  File  (AFRL  Model) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  AFRL  Model  control  By  Eigenstructure  Assignment  Technique 
%  Constraint  Function  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  [c,ceq]=nonlcon_e_val(x,A,B,C,max_k,lamd,vd,Wn) 


for  kk=l:9 

lamd(2*kk-l)=x(kk,l)+sqrt(Wn(kk)A2-x(kk,l)A2)*i; 
lamd(2*kk)=conj(lamd(2*kk- 1)); 
end 

%  Checking  the  result  of  assigning  the  eigenvalue. 
P=eye(20);  %eigenvector  penalty  matrix 


for  j  j = 1 : 20 

N=[lamd(jj)*eye(20)-A,  -B  ,zeros(20); 
zeros(9,20+9)  ,B'; 

P  ,zeros(20,9),  (lamd(jj)*eye(20)-A)']; 
VWN(:  jj)=inv(N)*[zeros(29,l);P*vd(:,jj)]; 
end 


for  jj=l:2:20 

VWN  ( :  ,jj  )=real(  VWN  ( :  ,jj)); 

VWN  ( :  ,jj + 1  )=imag(  VWN  ( :  ,jj + 1 )) ; 
end 

Kbar=-VWN(20+l:20+9,:)*pinv(eye(20)*VWN(l:20,:));  %  required  gain  matrix 


K_max=max(max(abs(Kbar))) ; 

c= [x( : ,  1 ) ;  K_max-  max_k] ; 
ceq=[]; 
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Appendix  L:  Objective  Function  File  (AFRL  Model) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  AFRL  Model  control  By  Eigenstructure  Assignment  Technique 
%  Objective  Function  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  f=fun_e_val(x,A,B,C,max_k,lamd,vd,Wn) 
for  kk=l:9 

lamd(2*kk-l)=x(kk,l)+sqrt(Wn(kk)A2-x(kk,l)A2)*i; 
lamd(2*kk)=conj(lamd(2*kk- 1)); 
end 

%  Checking  the  result  of  assigning  the  eigenvalue. 

P=eye(20);  %eigenvector  penalty  matrix 


for  j  j = 1 : 20 

N=[lamd(jj)*eye(20)-A,  -B  ,zeros(20); 
zeros(9,20+9)  ,B'; 

P  ,zeros(20,9),  (lamd(jj)*eye(20)-A)']; 
VWN(:Jj)=inv(N)*[zeros(29,l);P*vd(:,jj)]; 
end 


for  jj=l:2:20 

VWN  ( :  ,j  j  )=real(  VWN  ( :  ,jj)); 

VWN  ( :  ,jj + 1  )=imag(  VWN  (:  ,jj + 1 )) ; 
end 

Kbar=-VWN(20+l:20+9,:)*pinv(eye(20)*VWN(l:20,:));  %  required  gain 


E=eig(A-B*Kbar);E=real(E);I=find(E>0);s=size(I); 

C  ON =cond(  A-  B  *  Kbar) ; 
if  CON<10A9  &  s(l)==0 

load  XX_modred  %  Contains  variable  for  the  White  Gaussian  noise  intensity 
Sw=diag(x);  %  White  Gaussian  noise  intensity 

%Closed  loop  root  mean  square(lqr) 

S_x=lyap( A-  B  *  Kbar,  B*S w*  B ') ;  %Variance  of  x  matrix 
S_y=C*S_x*C';  %Variance  of  y  matrix 

RMS=diag(sqrt(S_y));  %Closed  loop  Root  mean  square  of  x 

f=RMS'*RMS/10A4; 


else 

f=10A20; 

end 
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Appendix  M:  Constraint  Function  File  (AFRL  Model) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  AFRL  Model  control  By  Eigenstructure  Assignment  Technique 

%  Constraint  Function  file 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  [c,ceq]=nonlcon_e_val(x,A,B,C,max_k,lamd,vd,Wn) 


for  kk=l:9 

lamd(2*kk-l)=x(kk,l)+sqrt(Wn(kk)A2-x(kk,l)A2)*i; 
lamd(2*kk)=conj(lamd(2*kk- 1)); 
end 

%  Checking  the  result  of  assigning  the  eigenvalue. 
P=eye(20);  %eigenvector  penalty  matrix 


for  j  j = 1 : 20 

N=[lamd(jj)*eye(20)-A,  -B  ,zeros(20); 
zeros(9,20+9)  ,B'; 

P  ,zeros(20,9),  (lamd(jj)*eye(20)-A)']; 
VWN(:  jj)=inv(N)*[zeros(29,l);P*vd(:,jj)]; 
end 


for  jj=l:2:20 

VWN  ( :  ,jj  )=real(  VWN  ( :  ,jj)); 

VWN  ( :  ,j  j + 1  )=imag(  VWN  ( :  ,j  j + 1 )) ; 
end 

Kbar=-VWN(20+l:20+9,:)*pinv(eye(20)*VWN(l:20,:));  %  required  gain  matrix 


K_max=max(max(abs(Kbar))) ; 

c= [x( : ,  1 ) ;  K_max-  max_k] ; 
ceq=[]; 


90 


Bibliography 


[1]  Karl  N.  Schrader,  Robert  H.  Fetner,  Steven  F.  Griffin,  Richard  Scott  Erwin 
Boeing-SVS,  Inc.  Development  of  a  sparse-aperture  testbed  for  opto¬ 
mechanical  control  of  space-deployable  structures,  Air  Force  Research 
Laboratory. 

[2]  Katsuhiko  Ogata.  Modern  Control  Engineering  Fourth  Edition ,  Prentice  Hall, 
Inc,  2002. 

[3]  Karl  N.  Schrader,  Robert  H.  Fetner,  Jed  Donaldson,  Robert  J.  Fuentes,  Richard 
Scott  Erwin  .Boeing-SVS,  Inc.;.  Integrated  control  system  development  for 
phasing  and  vibration  suppression  for  a  sparse -array  telescope.  Air  Force 
Research  Laboratory. 

[4]  Huckabone,  Thomas  C.  Algorithm  for  Robust  Eigenstructure  Assignment  Using 
the  LQR.  MS  thesis,  Air  Force  Institute  of  Technology,  Wright  Patterson 
AFB,OH,December  1991 

[5]  Meirovitch,  Leonard.  Element  of  Vibration  analysis.  New  York:  Me  Graw  Hill 
Publishing  Co,  Inc,  1986. 

[6]  Robinson,  J.  D.  A  Linear  Quadratic  Regulator  Weight  Selection  Algorithm  for 
Robust  Pole  Assignment.  MS  thesis,  Air  Force  Institute  of  Technology,  Wright 
Patterson  AFB,OH,  1990 

[7]  Jeffrey  B.  Burl.  Linear  Optimal  Control ,  Addison  Wesley  Longman,  Inc, 1999. 

[8]  Arther  E.  Bryson,  Jr.  Dynamic  Optimization ,  Addison  Wesley  Longman, 

Inc,  1999. 

[9]  David  Poole.  Linear  Algebra  ( A  Modern  Introduction  ),  Thomson  Learning, 

Inc, 2003. 

[10]  Harris,  Adam  G.  Robust  Eigenstructure  Assignment  Using  Positive  Definite 
Output  Feedback  Control.  MS  thesis,  Air  Force  Institute  of  Technology,  Wright 
Patterson  AFB,  OH,  December  1999. 

[11]  Richard.  Cobb,  Lecture  notes.  Air  Force  Institute  of  Technology,  Wright 
Patterson  AFB, OH 

[12]  Matlab  Control  System  Toolbox  User  Manual.  Lyap.  Limitation. 

Math  works  .2001. 

[13]  M.K.  Powers,  E.K.  Hall,  R.S.  Erwin,  “ UltraLITE  Deployable  Optical  Telescope 
(DOT)  Technical  Requirements  Document  -  Release  5.7,”  Air  Force  Research 
Laboratory,  AFRL- VS -TR- 200 1-1078,  May  9,  2001. 


91 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  074-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  the  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 

Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  an  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY) 

09-03-2004 

2.  REPORT  TYPE 

Master’s  Thesis 

3.  DATES  COVERED  (From  -  To) 

Sep  2002  -  Mar  2004 

4.  TITLE  AND  SUBTITLE 


5a.  CONTRACT  NUMBER 


6. 


CONTROLLING  THE  PRIMARY  MIRROR  IN  A  SPACE-BASED 
TELESCOPE  UTILIZING  AN  EIGENSTRUCTURE 
ASSIGNMENT  TECHNIQUE 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


AUTHOR(S) 


5d.  PROJECT  NUMBER 


Kim  Kyongsoo,  Captain,  ROK  Army 


7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(S) 

Air  Force  Institute  of  Technology 

Graduate  School  of  Engineering  and  Management  (AFIT/EN) 
2950  Hobson  Way,  Building  641 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFIT/GAE/ENY/04-M10 


9. 


WPAFB  OH  45433-7765 

SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFRL/VSSV 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


Dr.  Scott  Erwin 
3550  Aberdeen  Ave.  SE 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


Kirtland  AFB,  NM  87117-5776 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  feasibility  of  using  a  sparse  array  of  mirrors  in  a  satellite  system  is  being  evaluated.  This  study  focuses  on  the  control  laws 
necessary  for  achieving  the  minimization  of  the  effect  of  vibration  which  happens  at  the  primary  mirror  of  the  satellite,  by  using  the  eigenstructure 
assignment  technique.  White  Gaussian  noise  was  assumed  as  the  external  input.  AFRL  has  developed  a  79  state  model  with  9  control  inputs  and  9 
outputs,  3  sparse  array  mirrors,  and  a  base.  To  extend  the  theory,  a  simplified  model  with  2  mirrors  and  a  base  which  has  12  states,  4  control  inputs 
and  8  outputs  was  developed.  The  system  is  not  completely  controllable  nor  observable  because  there  are  4  rigid  body  modes . 

The  control  law  starts  from  realizing  the  states  which  are  controllable  or  uncontrollable.  Then  only  for  the  controllable  states,  the 
eigenstructure  assignment  technique  is  applied.  The  closed-loop  eigenstructure  was  determined  by  analyzing  the  open-loop  modes.  To  find  the  best 
value  for  the  eigenstructure,  Newton's  line  search  method  was  applied. 

The  relative  motions  of  the  primary  and  secondary  mirrors  were  described  by  wave  reflection  dynamics .  Results  show  that  this  technique 
achieves  a  reasonable  amount  of  control  use  for  the  objective  and  relates  these  result  to  those  obtained  from  linear  quadratic  regulator  (LQR) 
technique. 

15.  SUBJECT  TERMS 

Eigenstructure  Assignment  Technique,  LQR,  Lyapnov  Function,  Riccati  Equation, 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 


17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER 
OF 

PAGES 

103 


19a.  NAME  OF  RESPONSIBLE  PERSON 

Richard  G.  Cobb,  Major  USAF  (ENY) 

19b.  TELEPHONE  NUMBER  (Include  area  code) 


937-255-3069',  e-mail:  richard.cobb@afit.edu 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-1 8 


U 


U 


u 


